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1 Introduction 



Solar prominences consist of relatively cool, dense plasma that is suspended in the 
solar corona at heights up to 100 Mm above the chromosphere. They are observed as 
"filaments" on the solar disk, where they are seen in absorption in strong spectral lines 
(such as Ha) and in the Extreme Ultraviolet (EUV) continuum. When described as 
"prominences" they are seen above the solar limb, where they appear as bright features 
against the dark background. In this review the terms "filament" and "prominence" 
will be used interchangeably. The existence of cool, dense plasma suspended in the hot 
corona has b een a mystery ever since the fi r st observations of filamen ts and prominences 
(reviewed in iHiravamaL 19851 ; IZirkerl . Il989l ; iTandberg-Hanssenl . 1995). The discovery by 
iBabcock and Babcockl ( 1955T ) that solar filaments lie between different polarities of the 
Sun's magnetic field, provided one piece of the puzzl e by identifying the primary sourc e 
of support and constraint for the prominence mass jKippenhahn and Schluterl . Il957l ). 
However, the magnetic structure of prominences is still not fully understood, with many 
observations and theoretical models differing on the exact nature of the magnetic field. 
As a result the physical processes governing the origin and subsequent behavior of the 
prominence plasma remain a lively topic of debate. 

In recent years, high-resolution observations of prominences from the ground and 
from space have greatly improved our knowledge of the structure and dynamics of 
prominence plasmas. There have also been renewed efforts to measure the vector mag- 
netic fields in and below solar p rominences using spectro-polarimetric methods (see 
iLopez Ariste and Aulanierl. 12007). These observations have inspired a wide variety of 
theoretical models describing different aspects of the physics of prominences. The pur- 
pose of this review is to report on the progress that has been made in recent years, since 
the comprehensive review bv l TandbL.-S n^l Both state of the art space and 

ground based observations, along with models in four different areas of prominence re- 
search will be described. The four areas are chosen to give a broad overview of solar 
filaments, from the smallest structures currently observed, to their properties found on 
their scale size and finally filament properties on a global scale. These specific topics 
are: 



1. Prominence magnetic structure : The first topic describes the basic structure 
and morphology of solar filaments, including new high resolution observations and 
magnetic field measurements obtained from ground and space. The relationship of 
filaments to their underlying magnetic field along with recently developed static 3D 
magnetic field models, many of which are derived from observations, are described. 

2. The dynamics of prominence plasmas (formation and flows): The second 
topic considers the underlying physics and complex range of models that may de- 
scribe the origin and behavior of the filaments mass. Models include those based 
around magnetic forces, such as injection or levitation mechanisms, to thermal 
pressure force models of evaporation and condensation. In particular, significant 
new work on thermal non-equilibrium models is discussed. 

3. MHD waves in prominences: The third topic considers the significant advances 
that have been made in the last 10 years in observing and modeling the small-scale 
oscillations observed in filaments and prominences. The effect of flows and damping 
mechanisms on linear-MHD waves is discussed. 

4. The formation and large-scale patterns of filament channels: The final 
topic considers the properties of filament channels and filaments in a global con- 
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text. In particular Ha observations and theoretical models for the formation of 
filaments and filament channels are discussed, along with schemes used to catego- 
rize filaments. 

Other aspects of prominence research, including result s from space-based spe ctroscopy 



of solar prominences, are considered in the review by ( Labrosse et all 120091 . hereafter 



Paper I). We focus on the properties of non-eruptive prominences. 



2 Magnetic Structure of Filaments and Prominences 

In this section we first describe recent observations relevant to prominence magnetic 
structure, and then describe related modeling. 



2.1 Observations 
2.1.1 Basic Properties 

Filaments are always located above Polarity Inversion Lines (PILs), i.e., lines on the 
photosphere where the radial component B r of the magnetic field changes sign. Fila- 
ments can be found above PILs inside activity nests consisting of multiple bipolar pairs 
of spots ("active regions filaments"), at the border of active regions ("intermediate fila- 
ments"), and on the quiet Sun ("quiescent filaments"), including the polar crown. This 
nomenclature refers to their location on the Sun relative to underlying magnetic fields 
and not with respect to their internal plasma motions. Intermediate and active region 
filaments are located at sunspot latitude belts, however quiescent filaments may exist 
over all latitudes on the Sun. On the quiet sun the magnetic fields are concentrated 
into discrete network elements that are well separated from each other, with much 
weaker fields in between. Therefore, a PIL on the quiet sun is actually a zone of mixed 
polarity, and the precise location of t he PIL can best be defined by using spatially 
smoothed magnetograms ( Jonei, [2004). An example of a quiescent filament seen in Ha 



through different instruments, along with the underlying photospheric magnetic field 
distribution can be seen in Figure [1] In the next paragraph we first consider the basic 
structure of filaments as seen on the solar disk. Following this, the varying structures 
that are found when prominences are viewed above the limb are illustrated. 

Filaments typically consists of 3 structural components: a spine, barbs, and two ex- 
treme ends. The spine runs horizontally along the top of the filament, although there 
may be sections along the filament where the spine is nearly invisible. The barbs pro- 
trude from the side of the filament (see Figure 1(b) bottom panel) and when observed 
closer to the limb, using standard instrumentation with moderate spatial resolution, the 
barbs are seen to extend down from the spine to the chromosphere below. The barbs, 
as well as the ends of the filament (also called "legs" ) may be a c ollection of thread s 
that appear to terminate at a single point or at multiple points ( Lin et all . l2008al ). 



When viewing a quiescent filament on the solar disk at high resolution, Ha observa- 
tions indicate that each of t hese three structural components consist of thin thread-lik e 
structures (see examples in lMalherbd . ll9^ : lMartml Il998l : IPecseli and Engvoldl . EoOoT l . 
In Figure [5] these thin threads may be seen, where the observations are at the limit 
of present day resolution. The threads are found to have widths of about 200 km 
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2004/1 0/06, 10:06-11:09 UT, Themli/MTR 



2004/10/06, 10:06-11:09 UT, Themis/UTR 




2004/10/08. 08:35 DOT 



2004/10/06, 11:12 UT, SOHO/MDI 




(b) 



(c) 



Fig. 1 Ho filament observed on 6 th October 2004 from iDudik et all ll2008fl . (a) Meudon 
spectroheliograph in Ha line core at 09.35 UT. The box corresponds to the THEMIS/MTR 
field of view, (b) THEMIS/MTR (10:06-11:09 UT) and DOT (08:35 UT) observations in Ha 
line center, (c) THEMIS/MTR (10:06-11:09 UT) and SoHO/MDI (11;12 UT) longitudinal 
magnetograms, both saturate at ±40 G. 



( Lin et all l2005al lbL l2008allb|) and are not necessary aligned with the structural com- 
ponent of the filament to which they belong. Individual threads have lifetimes of only 
a few minutes, but the filament as a whole can live for many days. These thin threads 
are thought to be aligned with the local magnetic field. iLin et al. (2005a) find that for 
quiescent Hcv filaments it is not possible to associate the ends of individual filament 
threads with bright points in G-band images (such bright points are proxies for small 
kilogauss flux elements in the photosphere). They conclude that filament threads are 
located in thin bundles of field lines that are longer than the observed threads , i.e., only 
a fracti on of each field line is filled with absorbing plasma. In a recent paper. ILin et all 
(2008a) argue that at high resolution, all filaments ranging from active region to qui- 
escent are made up of such thin threads. While short threads are commonly seen in 
quiescent filaments, active region filaments seem to be composed of relatively longer 
threads. These differences may be related to the average angle made by the magnetic 
field to the axis of the filament spine. 

In contrast to the view on the disk, when observed above the solar limb differ- 
ent prominences can have very different appearance as can be seen from the three 
examples in Figure [3] In some cases the prominence consists of a collection of nearly 
horizontal threads or elongated blobs (Figure [3|a)), similar to the features seen in 



Fig. 2 High-resolution Ha images obtained from t he Swedish Solar T elescope illustrating thin 
threads aligned with the local mag netic field (from lLin et al\ 11200531) 1 . 



filaments on the disk. High-resolution Ha and Ca II H observations of prominences 
have rec ently been obtained w ith the Solar Optical Telescope (SOT) on the Hinode 
Satellite. lOkamoto et al. | (|2007h observed horizontal threads in a prominence near an 



active region and studied the oscillato ry motions of these th reads. A movie of the evo- 
lution of these threads may be seen in lOkamoto et all 1 20071 ). These horizontal threads 



are most likely aligned with the local magnetic field. In other cases, th e prominence 
consis ts of a collection of quasi- vertical threads (Figure E^b) and (c)). iBerger et al\ 
( 200sj) observed a hedge-row prominence and found that the prominence sheet is struc- 
tured by both bright quasi-vertical threads and dark inclusions. The bright structures 
are down- flow streams with velocity of about 10 km s , and the dark inclusions 
are highly dynamic up- flows with velocity of about 20 km s -1 . It is unclear what 
drives these up-flows and down-flows. The down-flow velocities are much less than the 
free -fall speed, indicating that the plasma is somehow bein g supported against grav- 
ity ( Pecseli and Engvold . I2OO0I ; iMackav fc Galsgaardl . l200ll ). How these quasi- vertical 



threads relate to th e magnetic field in such hedge-row prominences is not well under- 
stood. For example. ISchmieder et al ] (|2009l l have shown that dopplershift line-of-sight 
velocities in such threads can be of the same magnitude as observed vertical velocities. 
Assuming field-aligned flows, the authors argue for the existance of an oblique mag- 
netic field in 3D, that is projected on the plane of the sky. In addition, it is unclear 
how the lower altitude quasi-vertical threads seen above the limb relate to the barb 
threads seen in filaments on the disk. A major challenge of prominence research is to 
reconcile the different structures observed in filaments and prominences when seen on 
the disk versus limb. 

Filaments an d pro minences can also be observed in the H I Lyman lines (e.g.. 



Korendvke_etj^ .l200 j;lHeinzel et all l2001bl : IPatsourakos and Viail2002 l: ISchmieder et all 
20071 : iGunar et all 120071 ) and in the He II 304 A resonance line l|Wanal. Il999f ). These 
lines likely originate in the prominence-corona transition region (PCTR) . While promi- 
nences simultaneously observed in He II 304 A and Ha must lie in the same mag- 
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netic configuration, t hey commonly have different morphologies I Wang et all 1 199. 



iHeinzel et all boOlal ; lAulanier and Schmiedeil |2002| ; iDudik et all l200St ). Differences 
may include the observed width and whether the spine is visible or not. Such differ- 
ences may be explained by different formation mechanisms of these lines (see sections 
8.1 and 8.2 in Paper I) and their optical depths. Another question that arises from con- 
sidering different lines in prominences, is whether lines formed from ions and neutrals 
should exhibit the same prominence structure and if so what role ion-neutral coupling 
plays in this?. Observations of prominences with the SUMER instrument on board the 
Solar and Heliospheric Observatory (SoHO) satellite have shown that the Lyman lines 
have asymmetric profiles, indicating that (1) there are multiple threads along the line 
of sight, and (2) the t hread s move relative to each other with velocities of order 10 
km s" 1 |Gunar et all \200d ). For more detailed discussion of these observations, see 
sections 8 and 9 of Paper I. 

It should be kept in mind that the "filament" plasma seen in absorption on the 
disk may not be exactly the same as the "prominence" plasma seen above the limb 
some days earlier (East limb) or later (West limb). One complication in comparing 
filaments and prominences is that the prominence fine structure changes continually 
with time. Also, active region filaments, those filaments lying within the centers of 
activity complexes, generally lie at such low heights (less than lOMm, the approximate 
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Fig. 4 Schematic of a filament channel with fibrils which lie (1) anti-parallel to one-anothcr 
on either side of the PIL and (2) nearly parallel to the path of the PIL. The anti-parallel 
alignment indicates that the magnetic field (arrows) is dominantly horizontal and points in 
the same direction on either side of the channel. 



height of the spicule forest) that they are difficult to see at the limb. Figure [3] (a) 
and (c) illustrate prominences which lie at active latitudes, but can be seen above the 
spicular forest. They may be seen, as they are not active region filaments but rather 
"intermediate filaments" which lie on the borders of active regions and can reach much 
higher heights. Some prominences can be clearly seen in He II 304 A above the limb, 
but may not be so easily visi ble on the disk because of te mperature and d e nsity effects 
(see section 8.2 of Paper I). ISchmieder et al\ (2003) and ISchwartz et all (|200fj ) show 
that cool plasma (~ 10 4 K) with very low density could be present in the vicinity of 
filaments. This plasma is not visible in Ha but is detectable in EUV due to absorption 
of UV line radiation by the Lyman continuum (see Paper I). This may explain why 
prominences do not have the same aspect as filaments when crossing the limb. 



2.1.2 Filament Channels 

Filaments are located in filament channels, described as regions in the chro mosphere 



surrounding a PIL wher e the chromospheric fibrils are aligned with the PIL ( Martrea . 
Il96d : rGaizauskasl . fl998). These fib rils are i nterpr eted as giving the direction of the 
magnetic field in the chromosphere. iFoukall (|l971al lbT) noted that the fibrils emanating 
from magnetic elements in the channel show a streaming pattern that is opposite 
on the two sides of the channel (see Figure [4j). Given the magnetic polarity of these 
elements, Foukal deduced that the horizontal component of magn etic field must point 
in th e same direction on the two sides of the channel (also see iMartin et all 1 19921 . 
1994). This magnetic field along the PIL is believed to extend some height into the 
corona, and the filament is embedded in this field. The existence of this horizontal 
field within fi laments was c onfirmed by direct measurements using the Zeeman and 



Hanle effects iHvd or. 



ll998bl ). lMartin et al. 



Il965l; iRustl. Il9671 ; iLerov et all Il983l ; iLerovl. Il989l ; IZirker et all 



(1993) introduced the concept of chirality of filament channels. 



They classified filament channels as either "dextral" or "sinistral" depending on the 
direction of the axial component of the field as seen by an observer standing on the 
positive-polarity side of the channel (Figure [5]). To determine the chirality of filament 
channels, which do not necessarily contain a filament, high resolution Ha images (to 



+ + + 

DEXTRAL 



+ + + 

SINISTRAL 



Fig. 5 The chirality of a filament channel is defined in terms of the direction of the magnetic 
field along the channel (denoted by arrows) when viewed by an observer on the positive polarity 
side of the channel. For a dextral (sinistral) channel, the magnetic field points to the right (left). 



resolve individual chromospheric fibrils) are required. In strong field regions usually 
only Ha images are required. In contrast, for weak field regions where fibril patterns 
are not strong, magnetograms may also be used to aid the determination of chirality, by 
using them to determine the polarity of the m agnetic elements from which the fibrils 
extend from or go into. iMartin et al\ (|l994T l showed that channels in the northern 
hemisphere are predominantly dextral, while those i n the south are predominantly 
sinistral (also see lLerov et qZ.lll983l : IZirker et a71 . ll997l ). The origin of this hemispheric 
pattern will be further discussed in Sections 15. 1.21 and 15. 51 

As magnetic fields in filament channels are believed to extend to higher heights 
and the filament embedded in this field. Filaments may also be classified as being 
of "dextral" and "sinistral type". The chirality of filaments may be deduced, either, 
indir ectly from that of the channel or directly from magnetic field measurements (also 
sec iLerov et all 1 19831 ). To date no simultaneous studies comparing the chirality in 
filaments determined both indirectly and directly has taken place. 

Spectral diagnostics of the magnetic field orientation in a prominence observed 
with SUMER can also be used. Lyman lines have reversed/non-reversed profiles when 
the magnetic field o f the prominence is perp endicular/parallel to th e line of sight 
( Heinzel et 0/1120051 ). This has been tested bv lSchmieder et 

~d\ l|2007l ) with a round 

shaped filament crossi ng the limb. In the EUV, filament channels are much b roader than 
Hcv fila ment channels ( Heinzel et a/.l . [2001al ; [Aulanier and Schmiederl . [2002l ). lHeinzel et al\ 
(2001a) explain this enhanced width by absorption of EUV line radiation by the H I 
Lyman cont inuum; lower emission in the c hannel ("emissivity blocking") is another 
mechanism. lAulanier and Schmiederl ( 20021 ) used their magnetic flux tube models to 
deduce that magnetic field lines with concavities may be present in EUV filament chan- 
nels where plasma is cool. The existence of concave up magnetic fie ld lines has been 
tested recently by polarimetr ic measurements, in filam ent channels l|Litesl . 120051 ) and 
close to the feet of filaments ( Lopez Ariste et all \200(j ) . 



2.1.3 Filament Barbs 

When observed from above (i.e., at disk center), the filament barbs are seen to protrude 
at an acute angle with respect to th e long axis of the filament, like ramps off an 
elevated highway I Martin et aZ.L [T992). Therefore, filaments can be classified as either 
right-bearing or left-bearing depending on the directions of the barbs as seen from 
above. Martin and collaborators found a one-to-one correlation between the right/left- 
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bearing structure of the barbs and the chirality of the filament channel: filaments in 
dextral channels have right-bearing barbs, and those in sinistral channels have left- 
bearing bar bs. An important goal of filament modeling is to explain the cause of this 
relationship. iBernasconi et al. 1 (|2005h have developed software for automated detection 
of filament barbs and the chirality of filament channels. 

When observed away from disk center, the barbs are seen to be inclined structures 
th at extend f rom th e filament spine to the chromosphere below (e.g., SST observations 
of iLin et al.. 2008a). The vertical extent of a barb is much larger than the gravita- 
tional scale height of the prominence plasma (about 200 km), so the barbs cannot be 
static field-aligned structures; the plasma must somehow be supported against gravity. 
Therefore, a key issue is the orientation of the magnetic field in the barbs. This issue 
will be further discussed in section [2~2l 

The relationship between the filament barbs a nd the underlying ph otospheric mag- 
netic field has been discussed by many authors. iMartin et al. 1 1994 ) argue that the 
ends of the barbs of filaments on the disk are con n ected to weak magnetic fields in 
between the network elements. IMartin and Echolsl (jl994T ) propose that the filament 
barbs are anchored in parasitic (minority) polarity elements, i.e., weak magnetic fields 
with opposite polarity compared to the domina nt (majority) po l arity elements on the 
side of the filament where the barb is located. lAulanier et all (|l998T ) show that the 
barbs move i n accordance with t he changes of para sitic polarities d uring one day of 
observations. IWanel l| 19991 . l200lh . Ichae et all (|2005h and ILin et ail d2005ah find that 
the ends of the barbs are located very close to small-scale PILs between majority and 
mino rity polarit ies on the side of the filament. 

IWand (|200lf ) and IWang and Muglachl (|2007h proposed that flux cancellation be- 
tween the parasitic polarity and the neighboring dominant polarity plays a key role 
in the formation of filament barbs. Flux cancellation refers to the disappearance of 
photospheric magnetic f lux by mutual inte r action of opposite polarity fields in the 
photosphere ( Livi et all . 1 19851 : IMartin et all ll985T ) . It occurs everywhere on the Sun, 
but is particularly common at large-scale PILs on the quiet Sun where opposite polarity 
elements intermix and a zone of mixed polarity is created. At the edge of such a mixed- 
polarity zone, parasitic flux elements are likely to cancel against dominant polarity 
elements. IWand (|l999l . | 200 lh argued that magnetic reconnection accompanying photo- 
spheric flux cancellation i s the dominant mechanism for injecting mass into quiescent 
prominences. In contrast. iMartens and Zwaa n| (|200lh propose that downward-inclined 
barbs, linked to parasitic polarity elements, arise as a result of failed cancellation, i.e., 
the motion of flux elements across the PIL without cancellation. Changes of minor po- 
larities in filament channels lead to strong changes of Ha filaments and prove directly 
the relationship between minor polarities and barbs/ends of filaments. The disappear- 
ance and reappearance of a part of a filame nt has been directly rel ated to emergence 
of magnetic flux , follow ed by canceling flux (|Schmieder et aZlfeoOrJ) . 

IZirker et all ( 1998ah observed counter-streaming of plasma in a quiescent filament 
(see section 4, Paper I). Using observations at three wavelengths in the Ha line, they 
found evidence that matter flows along the spine of the filament in both directions. 
Similar counter-streaming occurs within the barbs: in one wing of the spectral line 
they observe up-flows from the barbs into the spine, while in the other wing they see 
down-flows. It is unclear what drives these flows, but they appear to be parallel to the 
magnetic field. Therefore, the observed flow direction can be used to infer the direction 
of the magnetic field, which may he lp resolve the questi on whether the magnetic fields 
in barbs are inclined or horizontal. IZirker et all 1 1998al ) propose that a filament barb 
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consists of a set of closely spaced flow channels that are highly inclined with respect 
to the vertical direction, with some channels having up-flows and others having down- 
flow. The velocities involved are 5 to 10 km/s, much le ss than the free-fall velocity 
corresponding to the fila ment heigh t. Deng et al. d2002|) found count er-streaming in 



an active region filament. iLin et al\ (|2003h and Schmieder et all 1 20081 ) observed such 



flows in large polar crown filaments two days before filament eruptions, so the counter- 
streaming may be a signature of the slow rise of a filament before eruption. With 
new high-resolution observations such as those by Hinode, many more examples of 
counterstreaming are expected to be found. With this it will be determined whether 
counterstreaming is a property common to all filaments and prominences, or whether 
it is a signature of the early stages of the onset of eruption. 

To investigate the nature of magnetic fields in the phot osphere below and a round 
barbs, polarimetric measurements were recently made by lLopez Ariste et~al\ (|2006h 
using THEMIS/MTR (Telescope Heliographique pour l'Etude du Magnetisme et des 
Instability Solaires). In Figure [S] an illustration of the 3D field can be seen, where 
white/black denotes the positive/negative normal field component and the red arrows, 
the direction of the horizontal field. Through using filament channel chirality rules to 
resolve the 180° ambiguity, the authors deduce that the field is of inverse polarity at 
the level of the photosphere. This type of photospheric configuration, with a collection 
of high altitude magnetic dips crossing the local PIL, as hinted from Figure [6] may fit 
many different structures of the coronal field lying above. An alternative interpretation 
of the coronal structure associated with the measured photospheric magnetic field is a 
vertical current sheet with no connections across the PIL. However such a structure is 
not stable to any form of non-ideal MHD perturbation and is likely to b e short lived . 
Non- ideal MHD perturbations would first result in the tearing of the field (|Furth et all 
Il963h and then the formation of an X-point structure (when viewed in a 2D projection 
along the PIL). Subsequent reconnection would then produce connections across the 
PIL. Depending on the force balance of the field and the nature of the extended current 
systems formed in the relaxation process, either a coronal arcade (of normal polarity) 
or bald patch (of inverse polarity) may result. As the long term existence of such a 
current sheet is unlikely, and since the photospheric field is of inverse polarity, the 
authors deduce that bald patches, i.e. magnetic dips that touch the photosphere, are 
present. These magnetic dips or hammocks are deduced to support the plasma in the 
barbs. Such observations are a major step in understanding the magnetic field structure 
in filaments and barbs. However to fully understand the magnetic field in filaments, 
new high resolution polarimetric measurements in many filaments are required. 

2.1.4 Measurements of Prominence Magnetic Field 

While observations of prominence fine structure provide some clues to the magnetic 
topology, the most direct determination of the prominen ce magnetic field come s from 
the inve r sion o f spectro-polarimetric data ( s ee rev iews bv lPaletou and Aulanierl . [2003; 
IPaletoul . 120081 : lLopez Ariste and Aulanierl 120071 1. Only a few spectral lines in the 
optical spectrum of prominences are suitable for magnetic field measurements using 
Hanle-effect diagnostic (the Hanle effect is the change in the polarization state of res- 
onantly scattered radiation due to prese nce of a magnetic field). The Ha line is strong 
and has interesting polarization profiles I Lopez Ariste et a/1 120051 ; iBianda et oflboOS : 



iRamelli et aUl^OOrj ), but the line is generally optically thick and the polarization trans- 
fer in this line is not understood. The lines most often used for measuring prominence 
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magnetic field are the He I multiples at 5876 A in the visible (He D3 ) and at 10830 A 
in the near-infrared. For recent attempts in modeling these lines see Leger fc Paletoul 
(2009|). The spatial resolution of such measurements is typically a few arc seconds, 
i.e., they cannot yet resolve the fine-scale structures seen in high-resolution images of 
prominences such as those obtained in SOT. 

A comprehensive effort to measure prominence magnetic fields was conducted in 
the 1970's and early 1980's, using the facilities at Pic du Midi (France) and Sacramento 
Peak Observatory (USA). The r esults of this work are review ed by iLerovl (1989), and 
can be summarized as follows l|Paletou and Aulanierl. 120031 ). The magnetic field in 
quiescent prominences has a strength of 3-15 G. The field is mostly horizontal and 
makes an acute angle of a bout 40° with respect to the long axis of the prominence 
( Bommier and Lerovl . [l998l ). The field strength increases slightly with height, indicating 
the presence of dipped field lines. Most prominences have inverse polarity, i.e., the 
component of magnetic field perpendicular to the prominence has a direction opposite 
to that of the potential field. The component along the pr ominence axis obeys a large- 
scale pattern jHvderlll965l;lRustlll967i;lLerov et aUI 19831) . which is now known as the 
chirality pattern IjMartin et all . 1 1994: Zirk er et aUll997h and will be further discussed 
in Sections 15.1.21 and 15.51 In the following discussion unless otherwise stated, angles 
given for the magnetic field in prominences are measured as an acute angle relative to 
the long axis of the prominence. Many features may effect these measurements, such 
as the position of the filament on the disk and the spatial resolution of the images. For 
each set of observations the reader should refer back to the original papers for such 
information. 

iLin et al\ ( 19981 ) presented the first full-Stokes spectro-po l arime tric observations of 
a filament on the disk in He I 10830 A. iTruiillo Bueno et al\ ( 20021 ) demonstrated the 
importance of lower-level atomic polarization in such measurements (the lower level 
of the 10830 A line is the ground level of the triplet system of neutral helium). They 
showed that this atomic polarization creates linear polarization of radiation in forward 
scattering at disk center. T his provides a new tec hnique for measuring the azimuth of 
magnetic fields in filaments. ICollados et al] ( 2003h applied this method to observations 
of two filaments in the Northern hemisphere. They find that the magnetic field is nearly 
horizontal, and that the magnetic vector is rotated with respect to the filament axis 
by an angle in the range 15 to 25 degrees. 
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IPaletou et al\ (|200lh reported full-Stokes observations of a prominence in He D3. 
Besides linear polarization signals due to the Hanle effect, they also found circular 
polarization signals an d they derived mag netic field strengths of 30-45 G in an active 
part of the prominence. ICasini et al\ 1 20031 ) published the first map of the vector field in 



a prominence. They fou nd that the average magnetic field in prominences is consistent 
with earlier studies (see iLerovl 1 19891 ) : a mostly horizontal field with a strength of about 
20 G and with the magnetic vector pointing 20° to 30° off the prominence axis. In 
addition, the map shows magnetic fields significantly stronger than average (up to 80 G) 
in clearly organized patches wi t hin the prominence (also see Wieh r and Biandal . [2003 ; 
lLopez Ariste and CasinJ . I2OO2L 120031 : lL6pez Ariste and Aulanj^jrH^OOTr i . It is unclear 
whether these structures are related to the fine structures of quiescent prominences 
( e . g., vertical th r eads). 

Casini et al ] l|2005l ) discuss various aspects of Stokes profile inversion, including 



the problem of finite optical depth of the He D3 line in the prominence, and the 
effects of line-of-sight integration of the polarization signals from different magnetic 
configurations. They confirm the presence of magnetic fields significantly stronger than 
average in limited regions of the prominence. They also emphasize the importance of 
full Stokes inversion for a correct diagnostic of magnetic fields in prominences. For 
certain values of magnetic field strength and inclination of the field with respect to the 
vertical direction, the Van Vleck effect produces a 90° ambiguity in the position angle 
of the magnetic field as projected onto the plane of the sky (in addition to the well- 
known 180° ambiguity). This 90° ambiguity makes it difficult to distinguish between 
horizontal and vertical fields in certain cases, and including Stokes V in the analysis 
helps disentangle this ambiguity. However, the results from the earlier studies based 
only on linear polarization measurements are thought not to be strongly affected by 
this ambigui ty (Lopez Ariste. priv ate communication). 

Recently. IMerenda et ail (|2006l ) observed He 1 10830 A in a polar crown prominence 



above the limb, and found evidence for fields of about 30 G that are oriented only 25° 
from the vertical direction. The reason for the different orientation of the magnetic field 
in this prominence (compared to the more horizontal fields found in other studies) is 
unclear. One possibility is the 90° ambiguity due to the Van Vleck effect, but the 
authors claim their particular measurement is unambiguous. 

2.1.5 Magnetic Fields and Flows in the Photosphere Below the Prominence 

An indirect method for probing the prominence magnetic field is to measu re the vector 
field in the photosphere below the prominence. lLopez Ariste et al\ 1 20061 ) deduced the 



presence of magnetic dips in the photosphere below a filament after having resolved 
the 180-degree ambiguity in photospheric vector-field measurements using the chirality 
rules for filaments (see Figure [SJ. These so-called bald patches are sites on the PIL 
where the horizontal field perpendicular to the PIL has "inverse" polarity, i.e., the 
field points from negative to positive polarity, opposite to the d irection expected for 
a potential field (see also iLitea . 120051 ). lLopez Ariste et al\ l|2006l ) argue that the dips 
are consistent with the presence of a weakly twisted flux rope in the corona above the 
PIL. If this is the case, the filament is thought to be supported by similar dips in the 
helical field lines at larger heights. 

Flows in the photosphere below the prominence may play an important role in the 
evolution and stability of the magnetic field supporting the prominence. Observations of 
shear flows, highlighted by maps of the vorticity, have been found below filaments, but 
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mostly during their eruption I Balasubramaniam et all 1 19981 : [Keil et all . 1 1999h . Until 



recently, observations o f orga nised photospheric flows below quiescent filaments were 
very rare. iRondi et al] ((2007) measured meso- and super-granular flows in the vicinity 



and below a filament before and during its eruption. They found that the disappearance 
of the filament starts where both parasitic and normal magnetic polarities are swept 
by a continuously diverging flow corresponding to supergranular diffusion. They also 
observed the interaction of opposite polarity elements, which is another candidate for 
triggering the filament destabilization as it implies a reorganization of the overlying 
coro nal magnetic f ield. 

iRoudier et all 1 20081 ) studied large-scale horizontal flows in the photosphere below 
a filament, using series of full-disk Dopplergrams and magnetograms obtained with 
the Michelson Doppler Imager (MDI) on SoHO. In Dopplergrams, the supergranular 
pattern is tracked using local correlation tracking (LCT), while in magnetograms LCT 
is used to obtain horizontal velocities. The topology of the observed large-scale flow- 
field changes during the eruptive phase suggests a coupling between the surface flows 
and the coronal magnetic field. An unusually fast north-south flow stream of amplitude 
30 to 40 m s" 1 was found to be compatible with the rotation of a part of the filament. 
This behavior suggests that the filament footpoints have been moved by the surface 
flows. The influence of the flow stream was enhanced by differential rotation. The 
authors measured an increase of the zonal shear in the region below the starting point 
of the eruption, and a decrease after the eruption. These results indicate that, even if 
large-scale flows do not lead to large-scale changes of the photospheric magnetic field, 
they nevertheless might be able to perturb the coronal field to an extent that filament 
eruptions and coronal mass ejections can be triggered. 

2.1.6 Coronal Structures Above the Filament Channel 

Quiescent prom inences are located in cavities at the base of coronal helmet streamers 



( Engvoldl .ri989). ii white-light observations above the limb, the cavity is darker than 



its surroundings a nd the height and width of the cav ity is roughly twice the height 



of the prominence ( Saito and Tandbcri 2,-Hansserllr973h . Some streamers contain con- 
centric arch systems (|Saito anTlTvd*e rril9*68T ) , suggesting the pr esence of an arcade of 



coronal loops overlying the cavity. lLow and HundhauserJ (| 19951 ) argue that the cavity 
represents a magnetic flux rope. 

In EUV and X-ray observations of filament channels on the solar disk, the coro- 
nal loops overlying a prominence often cannot be clearly identified, perhaps because 
they are too faint compared with the lower-altitude structures. However, X-ray ar- 
cades can often be seen in act ive regions, and after eruptive events on the quiet Sun 
( Martin and McAllisteri . Il996h . The coronal loops in these arcades do not cross the 



PIL at right angles, but are skewed in the direction along the PIL ("sigmoids"). The 
loops can be right-bearing or left-bearing with respect to the direction along the PIL 
(as seen from above). The observations indicate that dextral (sinistral) channels have 
left (right)-skewed coronal arcades, opposite to the right (left)-bearing structure of 
the filament barbs jMartin and McAllisteri . 1 19961 : iMartinl . Il998l : iMcAllister et all 120021 : 
ISchmieder et all . 12004 ). Taking into account the polarity of the dominant flux on ei- 
ther side of the filament, this implies that the axial component of magnetic field in 
the arcade is the same as that in the filament channel. Therefore, this implies that 
the filament channel observed in the chromosphere is part of a larger structure that 
extends to heights well above the prominence. 
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2.2 Models of Prominence Magnetic Structure 



The above-mentioned observations show that prominences are embedded in magnetic 
fields that are highly non-potential. As the dominant component of the field lies along 
the filaments long axis, such non-potential fields exhibit strong magnetic shear. How- 
ever, the detailed structure of these fields and the associated electric currents are not 
well understood. In the following we review various models for prominence magnetic 
structure. To begin with we focus on models that describe the structure as it exists at 
one instant of time, and we in general will ignore for the moment the question of how 
the field came to be that way. In Section Owe will see that the mass in prominences 
is infact dynamic. However since the flow speed is generally much less than the Alfven 
speed and within filaments the plasma-/? is low, such motions are often expected to 
be along field lines and not to strongly effect the field. Therefore, under this approxi- 
mation, static modeling is useful to describe the large-scale magnetic field structure of 
the prominence. This field is expected to evolve on times-scales much longer than that 
of the observed flows. Issues related to the formation of the prominence magnetic field 
and its evolution will be addressed in Sect. [5] 

An important question in prominence modeling concern s the role of magnetic dips 
in the structure and dynamics of the prominence plasma. iKippenhahn and Schluten 
(|l957l ) showed how cool plasma can be supported against gravity by a magnetic field 
containing dips in the field lines, i.e., sites where the field lines are locally horizontal 
and curved upward. For example, dips can occur in a potential field when the underly- 
ing photospheric magnetic sources have a quadrupolar structure ( "quadrupolar dips"). 
However, the magnetic fields in filament channels are not potential, and contain sig- 
nificant electric currents. Dips in non-potential magnetic fields can occur for several 
reasons: 



The structure of the magnetic field is significantly affected by the weight of the 
prominence plasma, distorting th e field in such a way as to create magnetic dips. 
Two examples are the isotherm al Kippenhahn and Schluterl (| 19571 ) model and the 
non- isothermal Fkood fc Anzerl ljl99(J ) extension. In both the shape of the dipped 
field lines is determined by the balance between magnetic and gravitational forces. 
Many authors have constructed magneto-hydrostatic ( MHS) models of prom i nence 
threads based on the Kippenhahn-Schluter model (e.g., Malherbe and Priestl . Il983l : 
iHeinzel and Anzerl. l200ll : iHeinzel et oZI. 120051 : ILqw and Petriel. l2005h Ano ther ex- 
ample is the "normal polarity" flux rope model bv lLow an^iZhana |2002l ). In this 
model the sense of magnetic twist in the flux rope is opposite to that in the overly- 
ing coronal arcade, and the force of gravity plays a crucial role in the existence of 
dipped field lines. In Paper I, Section 9.2 the Kippenhahn-Schluter model is used 
to compute the hydrogen spectrum in non-LTE modeling. 

Another possibility is that the dips exist in the absence of prominence plasma, i.e., 
the dips are not caused by the weight of the prominence. According to such models, 
the filament plasma is likely to be located near the dips of the field lines (see, how- 
ever, Section[3]). If the plasma /? < 1, these currents must flow nearly parallel to the 
field lines, so that the magnetic field is approximately force- free (fi is the ratio of gas 
pressure and magnetic pressure). Such force-free fields can contain dips in the field 
lines even when the underlying photospheric field is dipolar. In fact, many authors 
have suggested that filaments are s upported by nearly force-free flux ropes that 
lie horizontally above the PIL (e.g., iKuperus and Raadul . Il974l ; IPneuma 

d Il983l ; 
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van BallcEOoiicn and Martens! 1 19891: iPriest et al 1 Il989l: Rust and Kumaij . Il994l : 
Aulanier and Demoulirlll998l : lchae et all l200ll : iGibson and FanL 120061 ). An alter- 



native force-free configuration with dips, the sheared arcade , has also been s t udied 
as a description of the filament channel magnetic st ructure IjAntiochos et al. I ll99"l 
foeVore and Antiochosl l200fjl : lAulanier efail 120061 ). 

In quiescent prominences with low field strengths (B < 10 G) we may have a 
mixture of these two cases: the existence of the dips is due to the presence of a large- 
scale flux rope or sheared arcade (as in case 2), but the magnetic field near the dips 
is significantly distorted by the weight of the prominence plasma (/3 > 1) . In this case 
there are significant perpendicular currents near the dips, and the flux rope or sheared 
arcade is not force free. 

In the remainder of this section we consider the second category of models which 
describe the 3D structure of the magnetic field in the prominence, under the approx- 
imation that the weight of the prominence plasma does not play an essential role in 
the large-scale structure of the magnetic field. These models differ in terms of the 
requirement of dips for the formation and persistence of filaments. While the models 
differ in this sense, one common feature is that the magnetic field has a strong ax- 
ial component along the PIL. Such a field may be described as strongly "sheared". 
Magnetic shear either in the form of an arcade or a flux rope may build up th rough a 
number of process. Examples include: s hearing motions dAntiochos et al. , 1994 ), emer- 



gence of the upper part of a f lux rope ([Manchester et ali i 2004 : ATchontis and Torokl , 



l2008l : IPeVore and Antiochosl. I2OO0I : iFanl 120091) or shearing motions followed by con 



vergence and cancellation o f pho tospheric flux ( van Ballegooiien and Martens! . Il990l : 



iLitvinenko and Wheatlandl [2005). The origin of the axial fields in filaments will be 
fully discussed in Section [5.11 where a wide range of models and mechanisms for pro- 
ducing magnetic shear ar e discussed. Return i ng to the topic of static configurations, 
Section [2.2.11 discusses the lMartin and Echols! dl994T l Empirical Wire Model, which sug- 
gests that dips are not essential. Following this we describe two models in which dips 
play increasingly important roles. These are the sheared arcade model (Section 12. 2. 2\ 
and finally the linear and non-linear flux rope models (Section I2.2.3fl . 



2.2.1 Empirical "Wire Model" 

Through following a time series of Ha obs ervations of a quies c ent fi lament, along with 
photospheric magnetic field distributions. iMartin and Echols! (| 19941 ). proposed an em- 
pirical "wire model" for the magnetic field geometry in a filament. A key assumption of 
their model is that the fine scale structure observed in Ha lies parallel to the magnetic 
field. In the "wire-model" the filament plasma is located on magnetic arches that are 
highly sheared in the direction along the PIL. Although they are highly sheared, the 
arches do not contain dips. This feature distinguishes this model from others discussed 
in later sections. Some field lines run along the entire length of the filament and out- 
line the filaments "spine". Other shorter ones run partially along the spine, but spread 
out from it and connect down to minority polarity elements on either side of the PIL. 
These shorter structures represent the filament barbs. As neither the spine or barb 
field lines contain dips, the model requires the existence of other (non-magnetic) forces 
to act parallel to the inclined field lines to supported the cool dense plasma against 
gravity. Without such forces, the cool plasma would fall down to the chromosphere in 
a matter of minutes, which is not observed. The model however does not address this 
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Fig. 7 The Martin and Echols ( 1994) model of filament magnetic structure. The four panels 
show a dextral filament as seen from different perspectives. The green f ield lin es show the 
filament spine, and the blue and red field show the barbs. From [Lin et al 



question. In Sectio n [3] possible candidates for such forces will be discussed. While the 
iMartin and Echolsl (1 19941 ) model does not consider the origin of the sheared arcade, 
such a structure may be formed by shearing motions of the photospheric footpoints or 
by the emergence of the upper part of a flux rope as will be discussed in Section 12.2.21 

Th e Martin & Echols model was further discussed and developed by iLin et all 
( 200881 ). who presented the diagram shown in Figure [7] The panels illustrate the mag- 
netic field geometry deduced from four high resolution Ha images of filaments taken 
from different angles and including different classifications (active region-top left, in- 
termediate - top right and two quiescent filaments - bottom). In each image the spine 
is represented by the green field lines and the barbs by blue and red field lines. Fi- 
nally the end of the filament is denoted by the dashed lines entering the network fields 
at either end. From high resolution Ha images the authors deduce that thread like 
structures are common to all classifications of filaments and that in some examples the 
barbs consist of long thin threads that are inclined with respect to the horizontal. The 
authors conclude that such thin threads as the basic building blocks of all filaments 
no matter their classification. Subsequently the only difference between classifications 
may be the relative importance of the spine, barbs and ends. 

A k ey assumption of the Martin & Echols model that is illustrated by ILin et all 
(2008a) is that as the barbs, which are inclined field lines, spread out from the main 
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Fig. 8 Sheared arcade model for magnetic field in a filament channel: (a) Sheared field lines 
(red) with overlying potential loops (green), (b) Distribution of normal polarity (red) and 
inverse polarity (dark blue) magnetic dips. Each field line is plotted up to a height of 300 km 
above the dip to simulate th e effec t of the pressure scale height of the prominence plasma. 
Adapted from I Aulanier et all J2002T ). 

body and only connect to parasitic polarities in the photosphere. If this is the only al- 
lowed case it follows that dextral channels produce filaments with right-bearing barbs, 
and sinistral channels produce left-bearing barbs, consistent with observations. How- 
ever, as the authors consider barbs as inclined field lines, the model does not explain 
■why these inclined field lines cannot extend upwards from network fields on either side 
and join into the spine, thus producing the opposite bearing orientation of barbs on 
either side. Therefore, the model only provides a partial explanation for the observed 
correlation between the right/left-bearing structure of the barbs and the chirality of 
the filament channel. Later when discussing flux rope models (Section I2.2.3[) we will 
see that these models can provide a possible explanation for this structural asymmetry 
when magnetic dips are considered. 



2.2.2 Sheared Arcade Model 

A "sheared magnetic arcade" may be formed for example by shearing motions of the 



phot o spheric footpoints loc alized around the PIL I Antiochos et aZ.l . ll994l : lDeVore and Antiochosl . 
l200d : I Aulanier et ad lioO^ i. by the emergence of the upper part of a flux rope or by 
large-scale vortical motions. Numerical studies of this filament-channel formation mech- 
anism began with a simple bipole embedded in a larger scale background dipolar field, 
subjected to strong footpoint motions, parallel to the PIL and confine d to a narrow 



zone on either side ( DeVore and Antiochosl . l2000l : IPeVore et all 120051 ). These opp 



o- 

sitely directed flows drag the innermost portion of the bipole into a zone of weaker 
overlying field, yielding elongated, low-lying field lines that bulge upward at their less 
constrained ends and hence become dipped. As this stressed system relaxes, a variety 
of field-line paths are formed: some dipped field lines lean across the PIL, producing 
inverse polarity compared with the initial bipole, while others reconnect in the corona, 
producing weakly twisted field lines. Consequently the resulting filament channel is 
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neither fully normal nor fully inverse polarity. Fig. [8] illustrates the resulting magnetic 
structure. 



2.2.3 Flux Rope Models 

In the remainder of this section we consider models in which the weight of the promi- 
nence plasma does not play an essential role in the large-scale structure of the mag- 
netic field. Such models may be split into two categori es : weakly twisted flux r opes and 
highly twisted flux ropes (e.g lRust and Kumar] |l994 ); lokamoto et all (120081 ) ). While 



two categories exist, we focus on models containing weakly twisted flux ropes with dips 
in the magnetic field lines; the prominence plasma is assumed to be located at these 
dips. Weakly twisted flux ropes are only discussed as when non-erupting filaments or 
prominences are observed, they do not appear to be highly twisted. This is commonly 
observed in Ha images and through flows which are mainly horizontal. In Section [5] 
when considering models for filament formation a wider range of flux rope models will 
be discussed including those with strongly twisted flux ropes. Weakly twisted flux ropes 
are similar to sheared arcades, in the sense that the magnetic field is dominated by 
the axial component. One constraint on such models is that the flux rope must be ap- 
proximately in force balance with its surroundings, i.e., the configuration must be close 
to magneto-static equilibrium, and this equilibrium must be stable. Therefore, a key 
feature of such models is the existence of a coronal arcade overlying the flux rope. The 
coronal arcade provides the magnetic tension forces necessary to hold down the flux 
rope in the low corona. Although as previously stated, weakly twisted flux rope models 
are similar to sheared arcades, there is one key difference. For flux ropes models the 
flux rope and arcade are independent flux systems with a separatrix surface between 
them. In contrast, for a sheared arcade only a single flux system exists. 

The sense of twist within the flux rope is generally assumed to be consistent with 
the direction of the magnetic field in the overlying coronal arcade. Since the promi- 
nence plasma is located in the dips of the helical windings, the flux rope model predicts 
that the magnetic field at the prominence has inverse polarity wit h respect to the sur - 
rounding photospheric fields, consistent with Hanle measurements (|Lerov et all Il984h . 
Furthermore, the models predict that filament channels with dextral orientation of 
the axial field have left-helical flux ropes and left-skewed coronal arcades, and sinis- 
tral channels have right-helical flux ropes and right-skewed arcades. This is consistent 
with the observation that dextral (sinistral) channels have left (right)- s kewed arcades 
( Martin and McAllister] . Il99d ; iMcAllister et all l2002i ; ISchmieder et all |2004| ). 



If the effects of gas pressure and gravity are ignored, the magnetic field B(r) must 
be more or less force-free with electric currents flowing parallel or anti-parallel to the 
field lines: 

VxB = qB. (1) 

Here a(r) is the so-called torsion parameter, which is constant along the field lines 
(B • Va = 0). The degree of twist of the flux rope is a matter of debate. In some 
models the magnetic field lines inside the flux rope are highly twisted, but in others 
the central part of the flux rope is only weakly twisted and the field lines inside this 
core are nearly parallel to the flux rope axis. The key difference between such models 
is the spatial distribution of a(r). In models with strong twist, a(r) has its peak on 
the axis of the flux rope, whereas in models with weak twist a(r) has a hollow core 
distribution with the peak value of a(r) occurring at the interface between the flux 
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rope and its local surroundings (e.g.. lBobra et aZll2008l : ISu et aZll2009h . The observed 
field-aligned flows may be difficult to explain with highly twisted flux ropes. 



2.2-4 Linear Magneto-hydrostatic Models 

Aulanier and collaborators have developed 3D magnetic models of filaments by ex- 
trapolating photosphcri c magnetograms into the corona, assuming either linear force 



free fiel ds (LFFF, Aulanier and Demoulin , ~998l) or linear magneto-hydrostatic fields 



(LMHS , I Aulanier et q/.l. Il998i . Il999l . I2OO0I : 



Aulanier and Demoulinl . 120031 : iDudik et all 



2008). Such models produce helical flux ropes overlying the PIL, and the filament 
plasma is assumed to be located at dips in the helical field lines. The flux rope may 
be perturbed by the presence of magnetic elements in the photosphere. In areas where 
parasitic polarity elements are located close to the flux rope, the dips extend away from 
the main body of the filament, c onsistent with the observation that filament barbs ar e 
located near parasitic polarities (|Martin and Echolsl . ll994l ; [Martinl . ll99l ; IWang| . [200ll '). 
A key feature of the flux rope model is that the dips that represent the main body 
of the filament are all of inverse polarity. Therefore only minority polarity elements 
may perturb them to produce extended dips away from the main flux rope that rep- 
resent barbs. Due to the direction of the axial field in the flux rope, combined with 
the inverse polarity nature, such a perturbation always results in a right /left-bearing 
struc ture of barbs in dextral/sinistral orientations of the axial field (jMartin et qZ.1 . 1 1992 . 
1 1994 ). Therefore the flux rope model, along with any other model that exhibits these 
same properties, provides a natural explanation for the observed orientation of barbs. 
In addition the flux rope model resolves the problem of the opposite orientation of 
filament fine structures and the overlying loops, which is a real problem for the Martin 
& Echol's model. 

The modeling by Aulanier and collaborators predicts that the lowermost dips form 
bald patches on a part of the PIL separating the parasitic polarity from the surrounding 
dominant flux. Therefore, the "ends" of the barbs are located at the PIL, and the 
field lines passing through the barb are not actually rooted in the parasitic polarity. 
The model agrees with recent observations of bald patches in the photos pheric vector 
field in areas where filament barbs are located (|L6pez Ariste et all 2006). This means 
that the en ds of the filam ent barbs wou l d not b e rooted directly in the photosphere as 
proposed bv lWand (|200lh and lLin~ "oZl l|2005ah . The observations of magnetic dips are 



in contradiction with the idea of barbs rooted in parasitic polarities ( Lin et aU l2005a) 
and with the interpretation of th e observed counter-s treaming motions as plasma flows 
along fixed, inclined field lines ( Zirker et all fl99 8a') . More polarimetric observations 



should be done jointly with Hct high spatial resolution observations of fine structures in 
quiescent filaments to distinguish between the different magnetic models for the barbs. 
For active filaments and dynamic fibrils the situati on may be differ ent, and the fine 
structures probably end in the dominant polarities ( Lin et all . l2008bh . 



IDudik et all l|2008u used THEMIS /MTR observations of the photospheric vector 
field to construct a LMHS model of a filament. Figure [5] shows selected magnetic field 
lines and a comparison of the modeled dips with Hct observations. The authors find 
significant departures from translational invariance along the PIL. They show that the 
flux rope is split into two parts: one part is rooted in a network element (P2) on the 
positive polarity side of the filament, while the other part is present along the full 
length of the computational domain. This shows that flux ropes can have a complex 
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Fig. 9 Linear magneto- hydrostatic model of a filament observed with THEMIS (Dudik et al, 
2008). The left panels show selected magnetic field lines as seen (a) from the side of the 
computational domain, (b) by an observer on Earth. The yellow field lines are rooted within 
the domain, while the red and blue lines pass through the left and right boundaries. The 
right panels show the locations of dips in the field lines as seen from the same perspectives. 
Each field line is plotted up to a height of 300 km above the dip to simulate the effect of the 
pressure scale height of the prominence plasma. All panels show the vertical magnetic field 
on the photosphere, either as a greyscale image (upper panels) or as blue and purple contours 
(lower panels). Panels (b) and (d) also show the co-aligned He* image from THEMIS. 



magnetic structure, and may consist of multiple strands that are twisted together. The 
axial magnetic flux may vary along the length of the filament. 



2.2.5 Non-Linear Force-Free Field Models 



Several authors have developed non-linear force-free field (NLFFF) models that include 
the local environment of the flux rope. This is important because the overlying coronal 
arcade plays a key role in the equilibrium and stability of the flux rope. For example, 
IRegnier and Priest! (|2007h construct NLFFF models by "extrapolating" observed pho- 



tospheric vector fields into the corona (also see Reenier et al\ . 12002 : IRegnier and Amaril . 
120041 : IRegnier and Canfieldl . 120061 : ICanou et all l200fj| ). They consider two active re- 
gions: one a decaying active region with strong currents (AR 8151), the other a newly 
emerged active region with weak currents (AR 8210). For the old decaying active re- 
gion the connectivity and geometry of the NLFFF model include strong twist and 
strong shear, and are very different from the potential field model. The twisted flux 
bundles store magnetic energy and helicity high in the corona (about 50 Mm). The 
newly emerged region has a complex topology and the departure from a potential field 
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is small, but the excess magnetic energy is stored in the low corona and is enough to 
trigger powerful flares. 

An alternative method for constructing NLFFF models of filaments was developed 
bv lvan Ballegooiienl 1 20041 ). The technique involves inserting a flux rope into a potential 
field based on an observed photospheric magnetogram, and then evolving the field in 
time to an equilibrium state using magneto-frictional relaxation. The advantage of 
this method is that the axial flux <S> a xi an< A poloidal flux F po [ of the flux rope can be 
specified, and these parameters can be varied until a good fit to the observations is 
obtained. If the axial and poloidal fluxes are not too large, the relaxation results in a 
NLFFF describing both the flux rope and its surroundings. The method was applied 
to a filament observed at the Swedish Vacuum Solar Telescope (SVST) in June 1998. 
This U-shaped filament had a prominent barb that exhibited interesting fine structure 
and internal motions. It was found that the dipped field lines in the NLFFF model 
reproduce some of the observed features of the filament barb, but not all. The barb is 
attributed to the presence of a magnetic element directly below the filament, so it is 
di fficult to decide whe ther this element is "dominant" or "parasitic." 

iBobra et al. (2008) developed NLFFF models of two active regions using the above 
flux-rope insertion method. The models are constrained by Ha observations of the 
filaments and by TRACE observations of coronal loops in the overlying corona. The 
best fit to the observations is obtained for a model in which the flux rope is only 
weakly twisted; the peak value of a occurs at the boundary between the flux rope and 
its surroundings, not at the flux rope axis. Bobra et al. find that the axial fluxes of 
the flux ropes are close to the upper limit for stability of the force-free equilibrium; 
however, no major eruptions occurred in these active regions either several days before 
or after these observations. This suggests that active regions can release magnetic free 
energy gradually, and that the build-up of fr ee energy does n ot necessarily lead to large 
flares or coronal mass ejections. In contrast. ISu et all 1 2009T ) modeled a different active 
region and found that the flux rope present in that region has an axial flux that is 
well below the threshold for eruption. Clearly, observations of filaments can provide 
important constraints on the structure and stability of magnetic fields in active regions. 
For filaments between active regions and on the polar crown it is not known how close 
they are to the limit of stability. 

In more recent modeling of NLFFFs, iMackav et all (|2009h consider how the struc- 
ture of dips in a filament may vary as a single bipole polarity is advected towards 
the main body of the filament. The authors demonstrate how the bipole may result in 
the break up of the filament and a significant change in vertical extent of the dips. In 
contrast to linear force-free models they show that through considering a time-series 
of evolved non-linear force-free fields that barbs may form when either the minority or 
dominant polarity is advected towards the PIL. 

In the previous discussion the observed structure and morphology of solar filaments 
has been discussed along with recent developments in constructing static 3D models 
of their magnetic structure. In the next section the review considers the origin of the 
mass within filaments, along with the properties of observed flows. 



3 Plasma Structure and Dynamics 



Although the Sun's magnetic field clearly is the primary source of support and con- 
straint for the prominence mass, the underlying physics governing the origin and sub- 
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sequent behavior of this material remains a lively topic of debate. The observed char- 
acteristics of prominence plasmas are discussed in depth in Paper I (sections 3 and 
4) with a brief summary given below. The key observational constraints that must be 
satisfied by any successful plasma formation model are as follows: 

1. Active region prominences tend to be shorter (of order 10 Mm), short-lived (minutes 
to hours), and lower (<10 Mm), while those formed outside active regions can be 
hundreds of Mm long, persist through multiple solar rotations, and extend as high 
as 100 Mm. 

2. Prominence plasma is predominantly dynamic, exhibiting horizontal and /or verti 
cal motions of order 10-70 km s _1 feubota and Uesugil. Il986l : Schmicder et al 



1991; Zirker et al. 



2007; Bcrgcr et al 



1998al; IKucera et all 120031 ; iLin et all 120031 ; lOkamoto et 
20081 ; IChae et all I2008T ) . Q uiescent filaments can go through 



phases of enhanced i nterna l motion and activity ( Martini . 1 19731 : iTandberg-Hanssenl , 
1 19951; IKucera et all , 2006), altering the amount of cool material in the filament 



( Kilper et ai[ \2009) 



19381: Schri 


vei. 2C 


01), but 


(Chae et al. 


, 2000; 


Liu et all 



but also i s observed flowing up from the chromosphere 



Prominences appear to be composed of multiple knots and threads w ith widths 
down to the resolution limit (^100 km) and typical lengths of 3 - 20 Mm jLin et all 



2003l.l2005a! : lHeinzel and Anzerl.l20oi : lGunar et aZ.Ll2007l.l20oi ; lBerger et all\2004 ). 

5. Barbs — lateral extensions from the prominence spine toward patches of parasitic 
polarity with vertical extents much greater than the local gravitational scale height 
- apparently are c onnected to the chromosphere Sig^ EH; jg, 
1 19891 : (Martini . ll99Sl : IChae et adl2005h . 

As is discussed below, no single model can explain the entire complex range of charac- 
teristics at present. 

It was recognized over 30 years ago that the large prominence mass must come 
from th e chromosphere , because there is not enough plasma residing in the a mbient 
corona l|Pikel'nerl . Il97ll ; ISaito and Tandberg-Hanssenl Il973l ; IZirker et aH ll994T ) . Suffi- 
cient mass to explain prominences must be extracted from the chromosphere, either 
through magnetic forces, which inject or lift cool plasma directly into the corona, or 
through thermal pressure forces, which evaporate heated plasma that subsequently 
condenses into prominence knots or threads. Here we summarize the current status 
of these models for the prominence plasma, and discuss one well-studied example of 
an evaporation-based mechanism — the thermal non-equilibrium model — in greater 
detail. 



3.1 Injection Models 

According to the injection models, cool plasma is forced upward in filament-channel flux 
tubes with sufficient force to reach the observed heights of prominence plasma. Most 
injection models (see Fig. 1 10 p generally invoke reconnection low in the solar atmosphere 
as the pivotal cause of ma ss at prominence temperatures being expelled into the corona 
(|Wang| Il999l : I Chad . l200lh . This conjecture is largely motivated by the well-observed, 
but poo rly understood, connection between flux cancellation and filament channel for - 
mation l|van Ballegooiien and Martend , ll989l ; lMartlrlll998l ; IWang and Muglachl . 120071 '). 
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Fig. 10 Illustration of injection of prominence plasma (purple) by reconnection between 
minority-polarity bipole (blue lines) and preexisting filament-channel field (black lines). The 
polarity inversion line is dashed. 



Because injection necessarily produces unidirectional flows in each flux tube, counter- 
streaming could arise either through reconnection occurring at the bases of different 
flux tubes on opposite sides of the PIL or through sequential reconnection at each foot- 
point. Some injection models propose that the reconnection sites are at the PIL (e.g., 
Idia d . l2003l ) . while other s suggest tha t the jets originate at minority-polarity intrusions 
offset from the PIL (e.q.. lWanel.ll999l) . Many active-region prominences provide strong 
evidence for the type of jetting ex pected from reconnection events, and may well be ex- 



plained by injection 1 ChaeL 120031 ). However, active region prominences are short-lived, 



low-lying and highly variable, and are sufficiently distinct from quiet-Sun prominences 
that different physical mechanisms may be involved. Prominence barbs also might be 
explained by the injection of cool plasma by reconnec tion between the f i lament chan- 
nel field and sma ll bipoles emerging in the channel ijLitvinenkoL I2OO0I : IWanel . l200ll : 



IChae et all 120051 ) . Both large - and small-scale up-flows have been re corded in quiet- 
Sun hedge-row prominences (|Zirker et all 1 1994 iBerger et all [2OO8), although it is 



difficult to determine whether these flows originate at the photosphere due to spicules 
and other obscuring chromospheric activity. Recent high-resolutio n observations 



by 
il, 



SOT reveal strong jetting throughout the chromospheric network (|de Pontieu et al. 
2007); it is possible that this activity occurs within filament channels on field lines 
capable of hosting prominence material. 

On the other hand, it is unclear whether injection can account for all aspects of 
the observed dynamic evolution of quiet-Sun prominences, such as the frequency with 
which cool plasma appears suddenly in the corona and the predominantly horizontal, 
fine-scale counter-streaming flows. Furthermore, the extent to which reconnection can 
drive cool filament-channel material as high as 100 Mm in the corona, without also 
heating this plasma, has not been demonstrated in a realistic three-dimensional (3D) 
geomet ry with applicable energy sou rces and sinks. In simulations of Yohkoh X-ray 
jets bv lYokovama and Shibatal l|1995f ). reconnection produces hot jets directly, while 
the associated (but not co-spatial) cool jets are produced by compression of nearby 
open field. A similar mechanism might explain the association between flux cancel- 
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Fig. 11 Illustration of levitation of prominence plasma (purple) by U-loop emergence or post- 
reconnection relaxation associated with flux cancellation. 



lation and injection of cool plasma. Another unresolved issue is whether cancellation 
reconnection preferentially occurs in or below the chromosphere; if the interacting flux 
systems reconnect instead in the low corona (see below), then cool, dense plasma will 
not be lifted or injected directly. As is well-known, the fact that cancellation is observed 
in the photosphere or chromosphere does not nece s sarily indicate that the associated 
reconnection site is at the same level l Zwaanl . ll987l : lHarvev et all Il999h . The temper- 
ature minimum is thought by some to be the most favorable location for reconnection , 
because the Spitzer resistivity is highest there ( Roumeliotis and Moorel . 1 19931: Sturrock , 
199£ ; iLitvinenko Q or because of enhanced plasma turbulenc^bmi^l 

2002). In summary, the essential distinguishing features predicted by the injection 
model are that photospheric or chromospheric mass is injected with substantial speed 
at locations of flux cancellation/reconnection, either at or away from the PIL, and that 
this mass rises into the corona at or near its original cool temperature. 



3.2 Levitation Models 

Levitation models propose that cool plasma is lifted by rising magnetic fields at the 
PIL and transported transverse to the magnetic field (see Fig. 1 1 1 f) . In one class of lev- 
itation models, the filament-channel magnetic structure is a highly twisted flux rope 
that brings up cool plasma as the axis and lower portions of the rope emerge above 
the photosphere; thus, the prominence plasma should reside in the upward concavities 
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of the helical field ( Rust and Kumarl . ll994T l. However, such concave-upwar d formations 
have been observed only rarely in non-e r uptin g prominences (see, e.g.. iLitesl . liooj 
lLopez Ariste et~al. 1 120061 ; lOkamoto et~aH [2007). An alternat ive lifting mechan ism is 
the relaxation of magnetic fields during emergence of U loops (iDeng et aUbOQOh or af- 
ter reconnection associated with flux cancellation (|vanJ3aUeiro oijen^ L 1989^_ 



Priest et al lll996l:IOhver et aZ.l.ll999l: Litvinenko and Martini .! 1999: Gal seaard and Longbottoml . 



1999l ; lLitvinenko and Wheatlandll2005l ; IWelsch et aZlboOsf T Simulations of the reconnection- 
levitation mechanism have verified that rising field lines are indeed produced by recon- 
nection between bipolar systems, but significant work remains to prove that this mech- 
anism is responsible for the observed moti ons and properties of prominence plasma. De- 
tailed calculations of reconnecting bipoles l|Galsgaard and Longbottornl . ll99'9l ; lvon Rekowski and Hoodl . 
2008) predict significantly different plasma properties, depending on the dimensional- 
ity, the initial conditions, and the terms included in the energy equation. For example, 
IGalsgaard and Longbottoml ( 1999h found a 20% density contrast between the elevated 
material and the background corona, much l ess than the 2 orders of magni tude de- 
duced from prominence observations, whereas Ivon Rekowski and H ood (2008) predict 
that the uplifted plasma can exhibit X-ray-emitting coronal characteristics, filament- 
like photospheric properties, or a mixture of both (note, however, that the latter model 
is intended pri marily to explain coronal bri g ht po ints). The difference between the two 
results is that IGalsgaard and Longbottoml l|l999h due to computational reasons only 
considered a simple hydrostatic atmosphere with a uniform temperature profile and 
did not include a transition region. One common result is that most of the uplifted 
cool plasma drains along the rising field lines onto the chromosphere. Quiet-Sun promi- 
nences that resemble a central pillar from which cool material streams outward and 
downward are the best candidates for supporting this scenario. The extent to which 
this process can lift photospheric or chromospheric material as high as 100 Mm into 
the corona also has not been demonstrated. To date, none of these levitation models 
has allowed flux to retract beneath the photosphere during reconnection, thus preferen- 
tially favoring upward motions. Flux emergence simulations have shown that emerging 
U-loops and the associated photospheric plasma do not readily break through the pho- 
tosphere nor rise to coronal heights characteristic of quiet-Sun and intermed iate promi- 
nences, even without t he inhibiting presence of a preexi s ting coronal field jFanl . [2001 



Archontis et a . 



2004; Manchester et al., 2004; Magara, 2006; Galsgaard et al., 2007 



Magara et al. 2008; Arc hontis and Torokl [2008 V Although the U-loops of the emerg- 



ing tube do not rise to coronal heights, the process of flux emergence may still produce 
a coronal flux rope with dips. A flux rope may form through the reconfiguration of 
emerged sheared field lines that lie above the emerging tubes axis. The flux rop e may 



emerged sheared held lines that he above the emerging tubes axis, the nux rop 
be formed either by magne tic reconnection IjManchester et al. 1 120041 ; iMagaral . 



lArchontis and Torokl . 120081 1 or through helicity injection by torsional Alfven waves 



( Fanll2009h .' Once the tube is formed above the photosphere, it may rise to coronal 



heights, in principle dragging cool photospheric or chromospheric plasma with it. Such 
a process acting during flux emergence, is unlikely to explain quiescent or intermediate 
prominences but would be of importance to active region prominences. This will be fur- 
ther discussed in Section T5. 61 One effect not considered thus far in 3D flux-emergence 
calculations, the partial ionization of portions of the chromosphere and photosphere, 
might enable more rapid emergence as well as the emergence of more flux (and frozen-in 
plasma) into the corona I Leake and A rbor. 2006|). As discussed above for the injection 
models, it is unclear where reconnection between such flux systems would occur in 
the solar atmosphere; for example, the reconnection in the lGalsgaard and Longbottoml 
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Fig. 12 Illustration of an evaporation-condensation model: hot up-flows (red) driven by heat- 
ing localized above the footpoints evaporate chromospheric plasma that ultimately condenses 
in the corona as cool prominence material (purple). 



1999T ) study begins in the corona, at the apex of a null-null separator, so the bulk of 
the mass lifted at that location is coronal. An alternative mechanism proposed to lift 
and support prominence plasma in the c orona is upward-propagating, weakly damped 
MHD waves (|Pecseli and Engvoldl . I2OO0I ); recent o bservational evidence fo r significant 
Alfvenic perturbations in chromospheric structures l|de Pontieu et aZll2007l ) makes this 
an intriguing suggestion, but further quantitative work is needed to evaluate this mech- 
anism rigorously. In general, the key features of levitation models overlap significantly 
with those of injection models (see above); however, the levitated mass does not travel 
as far or as fast as injected mass would, and typically is predicted to be located above 
the PIL. 



3.3 Evaporation-Condensation Models 

Evaporation-condensation models are based on the fact that adding heat to a coronal 
loop increases the density of the corona while decreasing slightly the chromospheric 
mass (Fig. I12p . Different temporal combinations of heating and cool ing were studied, 
without success f ully reproducing the basic p roperties of a prominence (|Engvold and Jensenl . 
Il977l : lAnL Il985l : IPoland and Mariskal . [ 19861) . Subsequent research found that concen- 
trating coronal heating near the footpoints of a loop should produce a cool condensa- 
tion at or near the apex l|Serio et all Il98ll ; iMok et all Il990l ; lAntiochos and Klimchukl . 
ll99ll : lDahlburg et al.Ul998l ). Although recent analysis of TRACE observations of some 
coronal loo ps independently indicated that the heating is con centrated near the loop 
footpoints (|Aschwanden et all l200ll : IWinebarger et all |2002), the observational evi- 



dence for the spatial distribution of coronal heating is inconclusive at best. Confirma- 
tion of the basic principle behind this evaporation-condensation model for prominence 
plasma formation required the use of adaptive-mesh numerical simulations, to han- 
dle the rapid birth and subsequent evolution of a new, t hin transition region a t each 
interface between the loop and the cool condensation ( Antiochos et all 1 19991 ). The 
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Fig. 13 Prominence locations predicted by the thermal non-equilibrium model, assuming 
that the filament-channel magnetic structure is a sheared arcade. The white line on the grey 
photosphere is the PIL, and the boundaries of the surrounding filament channel are denoted by 
dashed lines. Only the blue field lines match the requirements for thermal non-equilibrium to 
produce dynamic and stationary condensations consistent with observed prominence structure. 



central concept of this model is that, if the heating scale A is small compared to the 
length of a coronal loop and localized near the chromospheric footpoints, then the 
plasma in the midsection of the tube, where the heating is negligible, must undergo a 
radiatively-driven thermal collapse to low temperatures. A runaway situation develops 
in the coronal plasma because the radiative losses increase with decreasing temperature 
T for T > 10 5 K as well as with the square of the density: once the local plasma has 
cooled to this critical transition-region temper ature, it must cool al l the way to chro- 
mospheric temperatures to regain equilibrium ( Hood &: Anzerj . Il98ct ) . The ratio of the 
heating scale to the loop length is a crucial factor because the total radiative losses from 
the loop increase linearly with length, but the thermal conduction and other energy 
transport or loss terms either decrease or remain constant with loop length. Clearly, the 
radiative losses will dominate for lengths above a thres hold value that is approximately 
an order of magnitude greater than the heating scale ( Mok et all . [l990;) . As discussed 
in section [2T2T the magnetic structur e containing the qu i et-Sun prominence material is 
most likely to be a sheared arcade llMartin and Echols. 1 199*3; iDeVore and Antiochosl . 
l2000h or weakly twisted flux rope (|Martens and Zwaanl . 1200 ll ; iBobra et all 120081 ) , in 
which many of the loops nearly aligned with the PIL are much longer than typical 
coronal loops. Therefore, for a given heating scale, condensations are more likely to 
form in these elongated loops than in shorter loops rooted outside the filament channel 
(see Fig. [I3). 
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A series of computational investigations of this evaporation-condensation process, 
denoted thermal non-equilibrium, has systematically explored the dynamics and ener- 
;etics of the plasma within individual elongated flux tubes heated near the footpoints 
Antiochos etail, Il999l . |2000| ; iKarpen et alt\200l\ , 12001 120051 , 120061 ; iKarpen fe Antiochosl . 



2008). While the nature of the footpoint heating is not well understood, this work has 
established constraints on the conditions favorable to condensation formation, the mag- 
netic structure of prominences, and the nature of the associated coronal heating. The 
key factors in determining the likelihood and behavior of condensations are the flux 
tube geometry and the localized heating properties. In long, low-lying flux tubes with 
shallow arches or dips, unequal quasi-steady heating yields a repetitive cycle of conden- 
sation formatio n, motion along the tube, and destruction by fall ing onto the nearest 
chromosphere ( Antiochos et all [2000; Kar pen et all l200ll . 120061 ). It is important to 



note that this process does not require the presence of dipp ed flux tubes, in con t rast t o 
the assumptions of many magnetic-structure models (e.g.. lAnzer and Heinzej . 120031 ). 
although shallow dips facilitate the collection and retention of cool plasma. On the other 
hand, if the host flux tube has a deep dip, as would occur in the outer portions of a 
highly twisted flux rope, thermal non-equilibrium produces condensations that rapidly 
fall to the lowest part of the dip and remain there, stationa ry but growing as long as 
the heating remains quasi-steady ( Karpen et fflZ.l . l2003l . 120051 ). Condensations also form 



when the energy input is impulsive i n nature and ran domly distributed in time, as in 
nanoflare models of coronal heating (|KlimchukL boOrt l. so long as the average interval 



between energy burst s is sh orter than the radiative loss time in the ambient corona 
(|Karpen fc Antiochosl . I2OO8I ) . The calculated condensation speeds, counter-streaming, 
lifetimes, and sizes are consistent with observations of many quiet-Sun prominences. 
However, this model does not provide a satisfactory explanation of active-region promi- 
nences, which are too short to support the thermal non-equilibrium process with typi- 
cal values of the heating scale, or of the vertical structure and dynamics of hedge-row 
prominences. Barbs could be consistent with thermal non-equilibri um if they are com- 



posed of vertically aligned dips in otherwise horizontal flux tubes ( Heinzel and Anzerl . 
bOOll ). as long as the dips are deep enough to trap the conden sed matter: for exam- 
ple, in photospheric "bald patches" near parasit i c polarity sites dAulanier et all [l998 ; 



lAulanier and Schmiederl 120021 ; Ivan Ballegooiienl . 120041 : lLopez Ariste et all \2006) 



This process could occur in any favorable magnetic structure, but thus far has 
been considered systematically only within the context of the sheared arcade model 
(section 12.20 . The color coding in Fig. 1131 illustrates the implications of the thermal 
non-equilibrium studies summarized above for the location and evolution of prominence 
material within a filament channel formed by the sheared arcade mechanism. Only the 
blue field lines match the requirements for thermal non-equilibrium to produce dynamic 
and stationary condensations consistent with observed prominence structure. Red field 
lines are too short; the helical black field line can host only a condensation sitting at 
the bottom of the dip; and green field lines are cap able of hosting sma ll , inte rmittent, 
short-lived condensations known as coronal rain ( Miiller et all 20031 . l2005h . Three- 



dimensional MHD simulations of both the magnetic and the plasma structure with 
localized heating, radiation, and thermal conduction are needed to determine the full 
scope and applicability of thermal non-equilibrium to solar prominence plasmas. For 
a discussion of energy balance considerations in solar prominences see section 10 of 
Paper I. Evaporation-condensation models, as represented by thermal non-equilibrium, 
uniquely predict both stationary and highly dynamic prominence threads that condense 
in situ in the corona and trace the supporting flux tubes; an increase in coronal density 
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precedes each condensation episode, while the collapse of the condensations reduces the 
ambient coronal density and generates waves and shocks. 

As discussed above, flows are commonly observed in solar filaments and promi- 
nences. In addition to this a wide variety of oscillations are observed. The properties of 
these oscillations, along with theoretical models used to describe them is now discussed. 

4 Magneto-hydrodynamic Waves in Solar Prominences 

4.1 Overview 

Theoretical studies of small-amplitude oscillations and MHD waves in the solar corona 
started long before clear observational evidence about the presence of these phenomena 
in the solar corona was available. The interest of their study lies in their potential 
relationship with the coronal heating problem and with the possibility to perform 
local seismology of the solar corona. Thanks to ground- and space-based observations, 
evidence about the presence of oscillations in coronal structures (loops, prominences, 
plumes, etc) is now widely available. 

Solar prominences are subject to various types of oscillatory motion (section 4, Pa- 
per I). In the early days of filament observations, all the data collected on oscillations 
were related with motions induced by disturbances coming from a nearby flare (e.g., 
iRamsev and Smith! Il966l ). These disturbances produce large-amplitude oscillations 
with velocity amplitudes of 20 km s _1 or higher. Observations of large-amplitude oscil- 
lations in filaments are rare, although in recent y ears with the help of new observational 
capabilities more detections have been repo rted (|Eto et qZ.Ll2002l ; IOkamoto et qZ.I . 12004 : 
Ijing et aLLl2003l . l200d : iGilbert efaH 120081 '). Also, theoretical modeling of such events is 
lacking. On the other hand, from observations performed with ground-based telescopes 
it is well known that most quiescent prominences and filaments display small-amplitude 
oscillations. The m otions are main ly detected through the periodic Doppler shifts of 
spectral lines (e.g.. rTsubakil . [l988h . and the observations have shown that (1) the os- 
cillations are of local nature; (2) simultaneously flowing and oscillating features are 
present; (3) the oscillations seem to be strongly damped in time, and (4) there are spa- 
tially coherent oscillations over large regions of prominences/filaments. A key aspect of 
prominence oscillations is the knowledge of its triggering mechanisms. While for large- 
amplitude oscillations they are well known, in the case of small-amplitude oscillations, 
the excitation process still remains a mystery, in spite of the available observational 
information. 

In this section, we will focus on small-amplitude oscillations in prominences, which 
can be interpreted as linear magneto-hydrodynamic (MHD) waves, and which con- 
stitute an important tool for understanding the physical properties (sections 2 and 
3 of Paper I) and the internal structure of prominences. Traditionally, the study of 
prominence oscillations has been based in the determination of the normal modes 
of oscillation of different equilibrium configurations, such as slabs or cylindrical flux 
tubes, without including any dissipative mechanism. Recently, the study of the damp- 
ing of prominence oscillations has been the subject of intense theoretical modeling by 
considering different dissipative mechanisms. Prominence seismology seeks to obtain 
information about prominence physical conditions from a comparison between obser- 
vations and theoretical models of oscillations. Attempts to determine the Alfven speed, 
the magnetic field, the shear angle, etc. have already been carried out. In the following 
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we first provide a summary of the observations of prominence oscillations, and then 
discuss various aspects of the theoretical modeling. 



4.2 Observational Background 
4-2.1 Ground-Based Observations 

Prominences above the limb exhibit oscillations with velocity amplitudes in the range 
2 to 10 km s" 1 . The oscillations appear to be local in nature, disturbing only some 
regions of prominences. Most information about this type of oscillation has been ob- 
tained from Doppler velocity data, although there are also detections coming from other 
spectral indicators such as line intensity and line width (section 2.1 and 3.4, Paper I). 
Only rarely are the oscillations detected in several of these spectral indicators at the 
same time and with the same period, which constitutes one of the puzzling features of 
prominence oscillations. Earlier, the first observational detections of oscillations using 
one-dimensional spectroscopic observations led to a classification in terms of short- and 
long-period oscillations. However, when the number of available observations increased, 
it was clear that a wide range of periods, between 1 and 90 minutes, were present. De- 
ta iled information on the oscillat ory p eriods detected in limb pro minences is available 
in lTsubakil (|l988h . lOliverl <|l999h and bliver and Ballesteil l|2002h . Although the clas- 



sification in terms of short- and long-period oscillations is still in use, it does not cast 
any light nor gives any help with regard to the nature, origin or exciter of oscillations. 

On the other hand, far more interesting r esults can be obtained from two-dimensional, 
high- resolution observations of prominences ( Molownv-Horas et qZ.LIl999l : lTerradas et all 



I2002T ) which allow the construction of maps of the Doppler velocity, wave period, damp- 
ing time and wave vector, and th e extraction of interes ting information about oscilla- 
tions in prominences. For instance. lTerradas et al\ [ 20021 ') reported the existence of large 



regions with periodic Doppler (line of sight) velocity oscillations having similar periods, 
around 75 minutes, noticing also that the oscillatory amplitude tends to decrease in 
time in such a way that the periodicity totally disappears after a few periods. Reliable 
values for the damping time, rn, have be en de rived from different D oppler velocity 
time series by iMolownv-Horas et al. I jl999h and iTerradas et 

Ml (|2002h . The values of 



Trj thus obtained are usually between 1 and 4 times the correspo nding period, and large 
regions of the prominence display similar damping times. Also, ITerradas et al\ (J2002) 
reported the presence, along two selected paths in the prominence region, of plane 
propagating waves as well as a standing wave. The plane waves propagate in opposite 
directions with wavelengths of 67,500 and 50,000 km and phase speeds of 15 km s _1 
and 12 km s _1 , respectively, while in the case of the standing wave the estimated 
wavelength is 44,000 km and a phase speed is 12 km s _1 . Furthermore, the analysis 
has identified the existence of a wave generating region, indicating that oscillations are 
locally excited. 

For the case of filaments on the disk, old t wo-dimensional observations of fi l ament 
illations |Thompson and Schmieden . Il99ll ; lYi and Engvoldl . Il99ll ; lYi et all \l99l\) 



oscil. 



Yi et al 



1991 



reveal that the Doppler signals of filaments form fibril-like structures, 
concluded that the fibrils form an angle of about 25° with the filament long axis, that 
individual fibrils or groups of fibrils may oscillate independently with their own periods, 
and that these fibril structures represent mass motions in magnetic flux tubes. Thanks 
to the improvement in spatial resolution provided by modern solar telescopes, high- 
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Table 1 Periods detected in prominence fine structure oscillations. 



Reference 


Period (min) 


Thompson and Schmicdcr (1991) 


4.4 


Yi, Engvold, and Keil (1991) 


5.3, 8.6, 15.8 


Lin (2004) 


26 


Lin et al. (2007) 


3 - 9 



resolution observations have given us key information about the internal structure and 
features of filaments, which is use ful for the study o f filament oscillations. For instance, 
observations of filaments made bv lLin et al\ (|2005al fbh suggest that the thickness of thin 
threads forming the filaments is < 250 km, that only a small portion of filament threads 
are filled with cool plasma, and that this absorbing pl asma is con tinuously flowing along 
the thread structure with velocities 15 ± 10 km s -1 . iLinl analyzed the Doppler 

signals obtained from two different regions of a polar crown filament, showing the 
presence of oscillation that are coherent over each observed region (size about 25 Mm) 
and are strongly damped after a few periods. Furthermore, in the middle part of the 
filament 49 Ha m oving features display periodic variations in Doppler velocity. Finally, 
iLin et al. | |2007l ) have found evidence for traveling waves in the threaded structure of 
a filament and, in some cases it seems that the propagating waves move in the same 
direction as mass flows. They also have determined the wavelength of a propagating 
wave, obtaining a value of around 3,000 km. In Table Q] the periods found from ground 
based observations are summarized. 



4-2.2 Space-Based Observations 



Prom inence oscillations have also been observed with space-based instruments. iBlanco et al\ 
(1999) used SUMER to study the behavior of different lines of Si IV and O IV in a limb 



prominence. The results show that a large amount of energy is contained in oscillations 
with periods between one and six minutes, corresponding to characteristic periods of 
the chromosphere and photosphere. Furthermore, for all the considered wavelengths, 
the energy versus the slit position varies in such a way that minima and maxima for 
each wavele ngth are coincident. L ater on, oscillations in an active-region filament were 
reported bv lReenier et al] 1 200 ll ). Using SUMER they observed the He I line, obtain- 



ing a time series of the filament with a duration of 7 h 30 min and with a temporal 
resolution of 30 s. A Doppler velocity time series was derived, and by Fourier analysis 
significant po wer was found at per iods between 5 and 65 min suggesting the presence 



of oscillations. iFoullon etaU l{2004 ) used SoHO/EIT 195 A with a time cadence of 12 
min to observe an EUV filament during its crossing of the solar disk. They reported 
intensity variations of long period, 8 - 27 h, with a dominant period of about 12.1 h, 
while the amplitude of the intensity variations reached 10% of the background intensity. 
Furthermore, the periodic intensity variations seem to be correlate d along the filament , 



and the most pronounced oscillations were detected during 6 davs. lPoueet et all (|2006 
used SoHO/CDS to observe two filaments, during different days, in the optically thick 
He I line (section 8.2, Paper I). The duration of the observations was 15 - 16 h with 
a time resolution of 20 s. Before performing the Fourier analysis of Doppler velocities, 
the mean velocity averaged over the width of the filament was computed and the ob- 
tained power spectrum points out the existence of a wide range of periods arriving up 
to 5 - 6 h. The results derived from these observations must be taken with care. First, 
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Table 2 Periods detected from space-based observations of prominence oscillations. 



Reference Period (min) 



Blanco et al. (1999) 


1 - 


6 


Regnier et al. (2001) 


5 - 


65 


Foullon et al. (2004) 


720 


Pouget et al. (2006) 


up 


to 360 


Okamoto et al. (2007) 


2 - 


4.5 


Berger et al. (2008) 


20 


- 40 


Ning et al. (2009) 


3.5 


- 6 



it should be noted that the He I line may be formed either by a 20,000K plasma or 
by scattering in a cooler 8,000K plasma. Therefore the emission may originate in the 
outer parts of the prominence, whereas the authors assume that the detected motions 
are propagated from the prominence core. Moreover, the averaging of Doppler signals 
coming from each pixel means that spatial information is lost and, furthermore, it is 
difficult to be sure that the same region of the filament is tracked when long-time 
observations are performed, and that there is no time evolution of the filament (bright- 
enings, fadings, small eruptions, etc). Recently. lOkamoto et al] ( 20071 ) using SOT and 
the Ca II H line have obtained about 1 h of continuous images of an active region. The 
Ca II H line movie shows continuous horizontal motions along the prominence. Some 
of the flows display constant velocities while others accelerate. Hinode movies also 
show that the threads suffer synchronous vertical oscillatory motions. Finally, recent 
high- resolution (0.2") observations of limb prominences made by SOT ( Berger et all 
2008) reveal very complex dynamics with vertical filamentary down-flows and vortices, 
as well as episodic, vertical up-flows. Furthermore, using horizontal time slices taken 
at different heights within the prominence, they suggest the presence of large-scale 
oscillations with periods between 20 and 40 min utes, lasting one o r two periods, and 
with a vertical phase speed of 10 km s _1 . Finally. iNing et "all ( 20091 ) have analysed the 
oscillatory behaviour of a quiescent prominence observed with Hinode. They find that 
prominence threads exhibit vertical and horizontal oscillatory motions. In some parts 
of the prominence, the threads seem to oscillate independently of each other, and the 
oscillations seem to be strongly damped. The periods reported are very short, with the 
dominant one appearing at 5 minutes. The range of periods determined from space 
based observations are summarised in Table [2] 

In summary, from the available observational information some characteristic fea- 
tures of prominence oscillations can be highlighted: The oscillations only affect parts 
of the prominence, and in the case of filaments are confined to the thread structure. 
The presence of simultaneously flowing and oscillating features is a common feature. 
The oscillations seem to be strongly damped in time. There are spatially coherent os- 
cillations over large regions of prominences and filaments. More detailed information 
about ob servation s of small amplitud e oscillations in filaments an d pro minences can be 
found in lEngvoldl |200ll . 12004 I2OO8I ), loiiver and Ballesteil (|2002r i, and lBaneriee et al\ 
J2007h. 



4.3 Theory of Small- Amplitude Oscillations in Prominences 



Theoretical interpretations of prominence oscillations are mostly based on linear MHD 
waves. An important ingredient in the modeling of MHD waves in prominences is the 
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Fig. 14 Sketch of the equilibrium configuration for a cylindrical model of a prominence thread. 
The grey zone represents the cold part of the loop, i.e., the prominence. The density in the 
prominence region is p p , in the evacuated (coronal) part of the loop, p e , and in the coronal 
environment, p c . The magnetic field is uniform and par allel to the z-axis , and the whole 
configuration is invariant in the (/^-direction. Adapted from Diaz et al. ( 2 0021) . 



chosen equilibrium configuration. In the past, slab models were used as a representa- 
tion of limb prominences. However, due to present day high-resolution observations, 
models now pay more attention to equilibria based on thin cartesian or cylindrical 
threads. These threads represent thin magnetic flux tubes, partially filled with cold 
plasma. Such structures seem to be the building blocks of filaments when observed on 
the disc. Furthermore, attempts have been also made to reproduce the observed damp- 
ing of prominence oscillations. The damping involves various thermal and non-thermal 
mechanisms. The thermal processes include damping due to optically thin radiative 
losses (section 8.1 and 10 of Paper I) and thermal conduction. The non-thermal pro- 
cesses include ion-neutral collisions in a partially ionized plasma (section 2.3, Paper 
I) and resonant absorption of wave energy in the inhomogeneous prominence plasma. 
In the following we first consider undamped, adiabatic MHD waves in prominence fine 
structures, and then discuss various damping mechanisms. Flows seem to be a charac- 
teristic feature of prominences (see section [3]), and the combined influence of flows and 
non-adiabatic waves on oscillations will be considered in section 14.3.51 



4-3.1 Adiabatic MHD Waves in Prominence Fine Structures 

Cylindrical geom etry seems to be the most suitable to model prominence threads. 
iDi'az et ali l|2002h considered a straight cylindrical flux tube with a cool region rep- 
resenting the prominence thread, which is confined by two symmetric hot regions 
(Fig. 1 14[) . In this case, the fundamental sausage mode (m — 0) and its harmonics 
are always leaky. However, for all other modes (m > 0), at least the fundamental 
mode lies below the cut-off frequency. Regarding the spatial structure of perturba- 
tions, in cylindrical geometry the modes are always confined in the dense part of the 
flux tube. Therefore, an oscillating cylindrical fibril is less likely to induce oscillations 
in i ts neighbor i ng fib rils, unless they are very close. 

iDfaz et~al\ 1 2005t) studied multithread systems in Cartesian geometry. In this case 



the equilibrium configuration consists of a collection of threads separated by a dis tance 
2c (see Fig. [To)) , and the individual threads are modeled as slabs ( Dfaz et aZ. , l200lh . Fol- 



lowing this approach, inhomogeneous filaments were constructed with different thread 
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Fig. 15 Sketch of the equilibrium configuration for a collection of fibrils. The grey zone 
represents the cold part of the loop, i.e., the prominence. The density in the prominence 
region is p p , in the evacuated (coronal) part of the loop, p c , and in the coronal environment, 
p c . The magnetic field is uniform and pa rallel to the 2- axi s, and the whole configuration is 
invariant in the y-direction. Adapted from Diaz et al. ( 2 0051 ) . 



density ratios representing the inhomogeneity in density of a real prominence, and 
the separation between threads has been chosen randomly within the realistic range. 
When the separation between threads is small, there is a strong interaction between 
them since the perturbation can easily overcome the separation. As a result, only one 
even non-leaky mode can cause all the threads to oscillate in phase. The dependence of 
the oscillation frequency on the thread separation c can be studied: in the limit c — s> oo, 
the structure and the frequencies of each thread oscillating alone are recovered, while 
for small values of c only the mode described before remains having a slightly smaller 
frequency than in the case of the single dominant thread mode. Therefore, for realistic 
values of the separation between threads, the multithread system would oscillate in 
phase, with similar amplitudes and the same frequency. These results, although de- 
rived in Cartesian geo metry, strongly ag ree with the observations of filament thread 
oscillations reported bv lLin et al . (2005a). However, in the case of cylindrical geometry, 
the spatial structure of perturbations indicates that the modes are always confined to 
the dense part of the flux tube. Therefore, an oscillating cylindrical fibril is less likely 
to induce oscillations in its neighboring fibrils, unless they are very close. 



4-3.2 Thermal Damping Mechanisms 



Soler et al. (2008a) studied the oscillatory modes of an equilibrium configuration rep- 



resenting a prominence thread (Fig. I16|) which is made of a homogeneous plasma layer 
with prominence conditions embedded in an unbounded corona. Thermal conduction 
parallel to the magnetic field, optically thin radiative losses, and heating have been con- 
sidered as non-ideal, damping mechanisms. Figure \T7\ shows the period P, the damping 
time tq, and their ratio, versus the longitudinal wavenumber, for the fundamental kink 
modes. Taking into account the results in the range of wavelengths typically observed 
in prominences, we observe that the internal slow mode (responsible for longitudinal 
motions) produces periods compatible with intermediate- and large-period oscillations, 
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Fig. 16 Sketch of the equilibrium configu ration for s t udying the damping of prominence 
oscillations in a single thread. Adapted from So ler et al ] ll2008ah 




to-'ior 1 jfl-*j 1 r*?ff- J 1<P to 1 

*„ fw-'J 





V. 




i 






p IIP 













10" 








Itp 
If 








X ,(,• 
10 s 




\ 

\ 




!£' 








fir* 


• 


■ 





m""nr" itr 1 to™ to-' to' to 1 

k (*>-'; 



to-"!/!™ ta- 1 10" in- 1 Hp 1tP 



Fig. 17 Period (left), damping time (center) and ratio of the damping time to the period 
(right) versus the longitudinal wavenumber for the fundamental kink oscillatory modes: internal 
slow (solid line), fast (dotted line) and external slow (dashed line). The shaded zones correspond 
to thos e wavelengths typically observed in prominence oscillations. Adapted from lSoler et ali 
(20083). 



whereas the fast mode (responsible for transverse motions) could be associated with 
short-period oscillations. On the other hand, the external slow mode mainly disturbs 
the surrounding corona, the amplitude of its motions within the prominence fibril being 
very small; hence it could be rather difficult to observe. Regarding the damping time, 
both internal and external slow modes are efficiency attenuated, with damping times 
of the order of their periods. However, again, the fast wave is much less attenuated 
since its damping time is between 2 and 6 orders of magnitude larger than its period. 

In order to assess the relative importance of each non-adiabatic damping mecha- 
nism, a comparison was made between the damping time obtained when considering 
all non-adiabatic effects (displayed in the middle column of Fig. I17|l and the results 
obtained when a specific mechanism is removed. The results of these computations sug- 
gest that (1) the internal slow mode is only affected by prominence-related mechanisms, 
radiative losses from the prominence plasma being responsible for the attenuation of 
this solution in the range of typically observed wavelengths; (2) the prominence thermal 
conduction is only efficient for very small wavelengths outside the observed range; (3) 
the fast mode is affected by both prominence and coronal mechanisms; (4) the damping 
of the external slow mode is entirely governed by coronal-related damping mechanisms, 
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mainly the coronal thermal conduction, which is the dominant mechanism in the range 
of typically observed wavelengths. 



4-3.3 Damping of MHD Waves in a Partially Ionized Prominence Plasma 

Prominences are partially ionized plasmas since hydrogen line s are observed (section 
1 an d 2.3 of Paper I). This fact was already considered by iMercier and Hevvaertsl 



(1977) when they studied, from a theo retical point o f view , the diffusion of neutral 



atoms due to gravity. A few years ago, iGilbert et all l|2002h studied the diffusion of 



neutral atoms in a partially ionized prominence plasma, con cluding that the los s time 
scale is much longer for hydrogen than for helium. Recently, IGilbert et al. I (|2007l l have 



investigated the temporal and spatial variations of the relative abundance of helium 
with respect to hydrogen in a sample of filaments. They have found that a majority 
of filaments show a deficit of helium in the top part while in the bottom part there is 
an excess. This seems to be due to the large loss time scale for neutral helium with 
respect to neutral hydrogen. 

The consideration of prominence plasmas as partially ionized is extremely impor- 
tant for the physics of prominences, and the effects on MHD waves in prominences need 
to be taken into account. In particular, the frictional damping of magneto-acoustic 
waves in a partially ionized plasma is much stronger than in a fully ionized plasma 
because the presence of neutrals causes the Joule dissipation to incre ase as a result of 
electron- neutral and ion-neutral collisions ( Khodachenko et all 12004 ) . A comparative 



study of the role of ion-neutral damping of MHD waves and their damping du e to vis- 
cosity and thermal conductivity was made bv lKhodachenko et all ||2004| . |2006| '). finding 
that collisional damping is dominant. 

The effects of partial ioniz ation on fast and slow waves in an unbounded prominence 
plasma have been studied by iForteza et al] ( 20071 ') . They have considered a partially 



ionized hydrogen plasma but have not included the effects of particle ionization and 
recombination, and they have assumed a strong thermal coupling between the species, 
so that electrons, ions, and neutrals have the same temperature (T c = T; = T n — T). 
The separate governing equations for the three species can be easily substituted by a set 
of one-fluid equations for the whole partially ionized plasma, where isotropic pressure 
has been assumed and gravity, viscosity, heat conduction, and non-adiabatic effects 
have been neglected. A uniform plasma with density po an d pressure po permeated by 
a magnetic field Bo = Bqx = const, has been considered, and linear magneto-acoustic 
waves have been studied. The results shown in Fig.[l8]suggest that ion-neutral collisions 
are more important for fast waves (for which t^/P is between 1 and 10 5 ) than for slow 
waves (for which rrj /P varies between 10 4 and 10 8 ). Therefore, one can conclude that 
the effects arising from the partial ionization of the plasma are unimportant for slow 
waves, while fast waves can be damped efficiently for moderate values of the ionization 
fraction, p, ~ 1, i.e., almost neutral plasmas. 

Another interesting feature is that when plotting the period or the damping time 
of the fast wave versus the wavenumber, at a certain critical wavenumber (fc c ) the fast 
wave disappears. The reason is that the critical wavenumber depends on the ionization 
fraction through r/Q (Cowling's magnetic diffusivity) . Hence, in a partially ionized 
plasma the fast mode only exists as a damped propagating wave for wavenumbers 
below the critical value, k c . For wavenumbers greater than this critical value we have 
a damped disturbance instead of a propagating wave. 
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Fig. 18 Period, damping time and Trj/P versus the ionization fraction for the fast wave (left) 
and the slow wave (right). The parameter values used are Tq = 8 000 K, po = 5 X 10~ 14 g/cm 3 
and Bq = 10 G, for which the Alfven speed is 126 km s _1 , the sound speed ranges from 10.5 
to 14.9 km s — 1 and the plasma j3 varies between 0.008 and 0.017. In addition, k x xo = 7r/2 
and k z xo = 0.1. The solid lines represent the numerical solutions of the dispersion relation 
while the triangles r epresent the results obtained with approximate expressions. Adapted from 
iForteza et~al\ J2007D . 



Apart from magneto-acoustic waves, the time damping of Al fven waves due to par- 
tial ionization effects has also been studied ( Forteza el q/.L lioO&t ) . Figure [19] shows the 
results obtained for the period, the damping time and the ratio of the damping time to 
the period. In this Figure the solution for a fully ionized plasma (jl = 0.5 ) with magnetic 
resistivity i Ferraro and Plumptonl . |l96ll : iKendall and Plumptonl . [l964h is also shown. 
The Alfven wave behavior is similar to that of the fast wave in the adiabatic partially 
ionized case. When ion-neutral collisions become the dominant mechanism, fast and 
Alfven waves have similar period and damping time and, as for the fast wave, the ratio 
of the damping time to the period decreases when going to almost neutral plasmas. A 
critical wavenumber (fcj) also appears for the Alfven waves. This quantity depends on 
the ionization fraction and on the propagation angle. Usually, k^ is larger than k c for 
fast waves, and both critical wavenumbers become equal for parallel propagation. 
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Fig. 19 Period, damping time and ratio of the damping time to the period for the Alfven 
wave in a partially ionized plasma with jl = 0.5 (dotted ) , jl = 0.6 (dashed), jl = 0.8 (solid) 
and jl = 0.99 (dash-dotted). Adapted from lForteza ei~al ] 1(20081) . 



4-3.4 Damping of Prominence Thread Oscillations by Resonant Absorption 



lArregui et al. (2008) studied the resonant absorption of transverse kink waves in fila- 



ment threads. They considered a straight, cylindrically symmetric flux tube of mean 
radius a, and neglected the effects of gravity. The inhomogeneous filament thread oc- 
cupies the full length of the tube, and is modeled as a density enhancement with a 
non-uniform radial distribution of density p(r) across the structure. The authors use 
the zero-/3 approximation, so that slow waves are absent, and they consider perturba- 
tions with azimuthal mode number m — 1, representing fast kink waves with transverse 
displacement of the tube. For fast kink waves to be damped by resonant absorption, the 
Alfven speed must vary across the structure. The authors assume that there is a tran- 
sition layer with thickness I between the interior of the filament and the surrounding 
corona. The ratio l/a provides a measure of the transverse inhomogeneity length scale, 
and can vary in between l/a — (homogeneous tube) and l/a — 2 (fully non- uniform 
tube). The global m — 1 kink mode is resonantly coupled to local Alfven waves. This 
coupling causes excitation of localized Alfvenic oscillations and damping of the global 
kink mode. 

Combining long- wavelength and thin boundary (l/a <C f) approxi mations, analytic 
expressions for the damping time over the period can be written as (e.g. .[flollweg and Yangl . 
Il988l : ISakurai et all Il99ll : boossens etaU Il992l . 1 19951 : iRuderman and Robert 1 120021 ) : 

l£ _ F £ £±i (2 ) 
P ~ b I c-1' 1 ' 

where c is the density contrast of the thread, and F is a num erical factor that dep ends 
on the details of the density variation in the transition layer. lArregui et al\ ( 20081 ) ob- 



tained numerical solutions of the MHD wave equations for the m = 1 modes. Figure I2TJ1 
shows that analytical and numerical solutions display the same qualitative behavior 
with density contrast c, wavelength A, and transverse inhomogeneity length scale l/a. 
For low values of the density contrast, the ratio r^/P rapidly decreases with increasing 
thread density (Fig. 120b ). Interestingly, for the larger contrasts typically found in fila- 
ment threads (c ~ 100), t^/P is nearly independent of c and A (Fig. 120b ). but rapidly 
decreases with increasing l/a (Fig. 120b ). 

The above results suggest that resonant absorption is a very efficient mechanism 
for the attenuation of fast kink waves in filament threads. It is not yet clear how this 
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Fig. 20 Ratio of damping time and wave period P for fast kink waves in filament threads 
with a = 100 km. In all plots solid curves correspond to analytic solutions given by equation 
J2} with F = 2/tt, and other curves correspond to numerical solutions, (a): Damping ratio 
as function of density contrast c with l/a = 0.2 and for two wavelengths, (b): Damping ratio 
as function of wavelength with l/a = 0.2 and for two density contrasts, (c): Damping ratio 
as function of transverse inhomogeneity length scale I /a for two combinations of wavelength 
and density contrast, (d): Percentage difference A of numerical and analytical values of the 
damping ratio, for A = 30a (dashed lines) and A = 200a (das h-dott e d line s), and for different 
values of density contrast (see labels). Adapted from lArregui ei~ai ] d2008f) . 



affects of overall dynamics of prominence plasmas. One possibility is that the resonant 
absorption causes significant heating of the prominence plasma, resulting in chromo- 
spheric evaporation (section 3). This may alter the density profi le of the prominenc e 
thread, which changes the Alfven wave resonance condition (e.q., lofman et aZlll998l ). 
Future modeling of the dynamics of prominence plasmas should include the effects of 
resonant absorption. 

4.3.5 The Effects of Material Flows 

As discussed in section [3] flows are a ubiquitous feature in prominences an d filaments, 
and are routinely observed in Ha, UV and EUV lines(see section 4, Paper D. ISoler et al\ 
(2008b) investigated the effects of both mass flow and non-adiabatic processes on the 
oscillations of an individual prominence thread. The thread is modeled as an infinite 
homogeneous cylinder with radius a, density p p , temperature T p ), and parallel flow 
velocity Up > 0. The cylinder is embedded in an unbounded homogeneous corona with 
density p c and temperature T c . The magnetic field is Bo everywhere and the total 
pressure is assumed to be continuous across the interface between the flux tube and 
the external medium. For simplicity, the hot coronal part of the magnetic tube that 
contains the thread is not taken into account, and gravity is neglected. The equilibrium 
configuration is shown in Figure [211 
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Fig. 21 Sketch of the equilibrium. Adapted from lSoler et al\ |2008b). 



Figure [22] shows the dependence of the period, damping time, and their ratio as a 
function of the flow veloci ty for the slow, fast an d thermal modes (for an explanation 
of the thermal mode, see ICarbonell et aZ|, boost) . The longitudinal wavenumber has 
been fixed to k z a = 10~ 2 , which corresponds to a value within the observed range 
of wavelengths. The flow velocities in the range 0-30 km s _1 are considered, which 
corresponds to the observed flow speeds in quiescent prominences. The anti-parallel 
slow wave becomes a backward wave for Up « 8.5 km s _1 , which corresponds to the 
non-adiabatic sound speed, causing the period of this solution to grow dramatically 
near such flow velocity. However, the period of both parallel and anti-parallel fast kink 
waves is only slightly modified with respect to the solution in the absence of flow, 
and the thermal wave now has a finite period, which is comparable to that of the 
parallel slow mode. The damping time of slow and thermal modes is independent of 
flow velocity, but the attenuation of the fast kink mode is affected by the flow. The 
larger the flow velocity, the more attenuated the parallel fast kink wave, whereas the 
opposite occurs for the anti-parallel solution. This be havior is due to wea k coupling of 
the fast modes to external slow modes (for details see ISoler et all l2008bh . 

Tcrradas et al modeled the tran sverse oscillations of flowing prominence 

threads as observed bv lOkamoto et ali (|2007h with SOT. The kink oscillations of a flux 
tube containing a flowing dense part, which represents the prominence material, were 
studied from both analytical and numerical points of view. The results determined 
that there is almost no difference between the oscillation periods when steady vs. 
flowing threads are considered. Also, the resulting period matches that of a kink mode. 
In addition, to obtain information about the Alfven speed in oscillating threads, a 
seismological analysis as described in 14.41 was performed. 



4.4 Prominence Seismology 

The main goal of prominence seismology is to infer the internal structure and properties 
of solar prominences from the study of their oscillations. One of the most searched for 
parameters is the prominence magnetic field strength and attempts to determine it have 



41 



a. 




11 , , , , , ! 10" I , , , , , I 10" I , , , , , 

5 10 15 20 25 30 5 10 15 20 25 30 5 10 15 20 25 30 

U r (km s-') U p (km ,-') U, (km s'') 




5 10 15 20 25 30 5 10 15 20 25 30 5 10 15 20 25 30 

U r (km s-') U, (km a'') U r (km s'') 



Fig. 22 Period (left), damping time (center), and ratio of the damping time to the period 
(right) versus the flow velocity for the fundamental oscillatory modes with k z a = 10~ 2 . Upper, 
mid, and lower panels correspond to the slow, fast kink, and thermal modes, respectively. 
Different line styles represent: parallel waves (solid line), ant i-parallel waves (das hed line), and 
solutions in the absence of flow (dotted line). Adapted from lSoler e t al. (20083). 



been made. For instance, using the significant periods obtained from the Fourier analy- 
sis of Doppler velocity time series belon ging to an active regio n filam ent, [Reamer et al\ 
(2001) applied the theoretical model of Joarder and Roberts! (jl993T i to determine the 
magnetic field and the angle of the magnetic field with the long axis of the observed 
filaments. They used the periods to identify the modes involved in the oscillations 
and after this identification they obtained 18° ± 2.5° for the angle, assuming a promi- 
nence temperature of 800 K and an analytic al relationship between the magnetic field 
strength and the density. IPoueet et all hood ) followed the same procedure in the case 
of filaments observed with CDS. They identified six modes (slow, fast and Alfven) in- 
volved in the oscillations and, using the same theoretical model, they obtained values 
between 19° and 35° for the angle and between 10 and 35 G for the magnetic field 
strength. All these determinations must be considered with care because the identifi- 
cation of the modes based on the periods is very uncertain; apart from the periods, the 
velocity polarization is of paramount importance for the mode identification in slab 

models. 

As discussed in section 14.3.11 iDiaz et~al\ d2002l ) considered a cylindrical model of a 
prominence thread. According to this model, the cool prominence thread has a radius 
b, and its length 2W is less than the length 2L of the flux tube on which the thread 
is located (see Fig. I14|l . The density in the prominence is p p , the density in the coro- 
nal part of the flux tube is p e , and t he density in the surrounding corona is p c . An 
important result of iDfaz et all (|2002l ) is that the dimensionless oscillation frequency 
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Fig. 23 Frequency of (a) and (c) the kink even modes and (b) and (d) the kink odd modes 
against p p /p c for the parameters pe/pc = 0.6 and b/L = 0.001. In (a) and (b) W/L = 0.1, 
while in (c) and (d) W/L = 0.2. The right axis provides with the period To obtained after 
assuming the background magnetic field strength, the coronal density and the half-length of 
field lines are Bo = 5 G, p c = 8.37 X 10 -13 kg/m 3 and L = 5 X 10 4 km . The horizontal dashed 
line in each plot gives the cut-off frequency. Adapted from lDfaz et all 1120021 ). 



ujL/cac depends mainly on the ratios pp/pc and W/L, and is insensitive to the ratios 
b/L and pe/pc- Here ca c is the coronal Alfven speed. Fig. 1231 shows the dimension- 
less frequency of some modes for fixed values of fibril length W/L, thickness b/L, and 
pe/pc- Using such diagrams, we can perform a seismological analysis if some of the 
involved quantities are provided by observations. The availability of this observational 
information would allow us to determine the numerical values of the rest of promi- 
nence physical parameters. These examples point out that in order to determine the 
magnetic field strength, and apart of the oscillatory period, other physical properties 
of the prominence, such as the density, must be known or assumed. 

In the theoretical model used bv lTerradas et al\ |200Sl ) (section l4~33j) . prominence 
and coronal Alfven speeds are related through the ratio between prominence and coro- 
nal densities. Therefore, once this ratio is fixed and the coronal Alfven speed is deter- 
mined from the dispersion relation for the kink mode, the Alfven speed in prominence 
threads can be determined. Figure [2H shows the dependence of the Alfven velocity in 
the thread as a function of the coronal Alfven velocity, for two different threads, when 
different density ratios and lengths of the magnetic field lines are considered. It can 
be seen that for a certain value of the coronal Alfven speed the prominence Alfven 
speed stabilizes. Using this analysis, lower bounds on the Alfven speeds in different 
prominence threads can be obtained. 



43 



Thread §2 




o \ 1 

1000 2000 3000 4000 

v Ac (km s~ } ) 

Fig. 24 Dependence of the Alfven veloc ity in the thread as a f unction of the coronal Alfven 
velocity for two of the threads studied by Okamoto et al. (2007]). In each panel, from bottom 
to top, the curves correspond to a length of magnetic field lines of 100,000 km, 150,000 km, 
200,000 km, and 250,000 km, respectively. Asterisks, diamonds, triangles, and squares corre- 
spond to density rati os of the thread to the coronal gas p p /p c — 5, 50, 100, 200. Adapted from 
iTerradas et al\ J200SI) . 



For the threads considered and taking the length of the magnetic field lines to 
be 100,000 km, the lower bounds of the Alfven speed are between 120 - 350 km s -1 . 
These values are consistent with strong magnetic fields (50 G) and large densities 
(1 — 8 x 10 11 cm -3 ), or weaker magnetic fields (10 G) and lower densities (0.4 — 3 x 10 10 
cm -3 ). These densities are within the ranges inferred from spectroscopic techniques 
(see sections 2.2 and 3.3 of Paper I). Further assumptions about the density, appropriate 
for the considered model, but not based upon simultaneous observations (which are 
not available), must be made for a numerical determination of the magnetic field. 

Finally, an important conclusion in lArreeui et cH\ l|2008h (section EXI)) is that 
in equation @ the damping rate becomes independent of density contrast for large 
values of this parameter. This fact has several seismological implications: First, the 
determination of the density contrast in the case of prominences is not as crucial 
as in other problems where the density contrast is low; second, by assuming a density 
contrast c of infinity in equation (0 an estimate of the transverse inhomogeneity length 
scale l/a can be determined. Using t^/P = 4, we find l/a ~ 0.15, i.e. high-density 
threads are compatible with thin inhomogeneous layers. Furthermore, knowing the 
period and the wavelength, the Alfven speed in threads can be also obtained. To make 
further progress and to obtain the magnetic field strength the prominence density is 
required. 



In summary, the determination of the magnetic field strength is not straightforward. 
Depending on the theoretical model assumed to interpret the oscillations, additional in- 
formation on numerical values of parameters such as prominence and coronal densities, 
geometrical dimensions, etc., need to be observationally acquired before the magnetic 
field strength can be determined. 

4-4-1 Concluding Remarks 

The study and understanding of small-amplitude prominence oscillations is a challeng- 
ing task from the observational and theoretical point of view. The number of different 
physical effects (non-adiabaticity, flows, partial ionization, etc.) involved suggests that 
they must be explored in a systematic manner in order to obtain a full understanding of 
the properties of these oscillations. From a theoretical point of view, the main question 
could be: What do we need in order to foster theoretical understanding of prominence 
oscillations? Of course, a simple answer would be better observations having high 
spatial and time resolution. However, this would not be enough because our current 
knowledge of the internal structure and physical properties of solar prominences is not 
complete. For example, the magnetic field structure which supports and shields the 
prominence is not fully understood. Without a clear understanding of the magnetic 
structure it becomes difficult to advance our knowledge of prominence oscillations. 
Furthermore, are the flows real or are there ionization or excitation fronts changing 
the physical conditions and mimicking flows? What is the mechanism responsible for 
prominence heating? Partial ionization seems to be a key feature in prominences, in- 
fluencing not only magneto-hydrodynamic waves but also the manner in which promi- 
nence plasma is tied to magnetic field lines and, right now, is starting to be considered 
in some theoretical studies. 

From the above considerations, one can easily conclude that from the theoretical 
and observational point of view we are far from a complete understanding of promi- 
nence oscillations, and this conclusion is correct. In the near future, the joint use of 
two-dimensional, high-resolution observations, modern data analysis techniques, and 
complex theoretical models, incorporating as much physics as possible and using state 
of the art numerical simulations, should help us to make further progress in this field. 



5 Formation and Large-Scale Patterns of Filament Channel and Filaments 

5.1 Global Patterns and Formation Locations 

Solar filaments (a.k.a. prominences) form over a wide range of latitudes on the Sun. 
Their locations spread everywhere, from the active belts to the polar crown. Poleward 
transport of magnetic flux across the solar surface during the solar cycle is ac companied 



by a poleward migration of the preferred locations of filament formation (|McIntosh 



2002: lMinaroviech et aZ.lll998l : lAmbroz and Schrolil20ollMouradian and Soru-Escautl . 



19941 ). Although filaments may form at many locations on the Sun, they always form 
above Polarity Inversion Lines (PILs), which divide regions of positive and negative 
flux in the photosphere. As discussed in Section 12.1.21 a necessary condition for the 
formation of a filament is the pres ence of a filament channel at the height of the chro- 
mosphere (e.g. , [Gaizauskasl . fl998 L Filament channels are the unique sites of filament 
formation. Channels are more fundamental than the filaments that form within them, 
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as not every channel contains a filament and a single channel may survive a succession 
of filament formations and eruptions. The structure of a filament channel is illustrated 
in Figure [3] A key feature of a filament channel, is that it is a region of dominant 
horizontal field where the field on either side points in the same direction. As a result, 
filament channels are interpret as locations of strong magnetic shear and highly non- 
potential magnetic fields. The channels are believed to be a low atmosphere signature 
of a larger horizontal non-potential field that extends into the coronal volume (see 
Section f2 . 1 . 2 1) . Presently it is unclear why channels and their non-potential fields build 
up along PILs, so understanding filament channel formation is key to understanding 
the evolution of magnetic fields on the Sun and their relationship to eruptive phenom- 
ena. By observing and interpreting their formation and evolution, we may examine 
directly the buildup to, and initiation of geoeffective space weather. In the discussion 
below we consider the global properties and formation locations as deduced from Ha 
observations. 



5.1.1 Global Classification Schemes 



While filaments form at many locations on the Sun, very few studies have considered 
the exact nature or history of the PILs above which they form. Those studies that 
have considered this, are mainly restricted to studying large-scale, stable filaments 
and neglect smaller unstable filaments forming in the centers of activity complexes. 
Understanding the type of magnetic environment in which filaments form is key to 
understanding the magnetic interactions required for their formation. 

Over the ye ars many cl assificatio n schem e s for filaments and prominences have 
been developed (d' Azambuia and d'Azambuial . 1 19481 ) . These schemes have referred to 
different features, including: the dynamics of the material in the prominence, whether 
it is stable or eruptive, or finally with re spect to the distribution o f the magnetic flux 
below the filament. For a discussion see iTandberg-Hanssenl (1995). In this paper we 
use a clas sification scheme of filaments in terms of their spatial location on the Sun 
( Engvoldl . ll998T ). In later sections when we discuss the possible mechanisms of fila- 
ment formation this classification scheme will prove useful in illustrating that different 
mechanisms may form different types of filaments. As given in the introduction the 
types of filament were described as: ARF (Active Region Filament), IF (Intermediate 
Filament), and QF (Quiescent Filament). An ARF is one which forms in the centers of 
active regions or activity complexes. In contrast, IFs form between active regions and 
decaying regions of unipolar plage. Finally QFs form in regions of weak background 
fields. Observations tend to show that IFs and QFs are larger, much more stable struc- 
tures with longer lifetimes (weeks to months) compared to ARFs, which are generally 
unstable with a lifetime of only a few h ours to days. 

While the classification scheme of lEngvoldl (|l998T l provides a useful distinction 
between filaments forming inside and outside active regions, to understand the role 
that magnetic fields play in the formation, structure and evolution of filaments it is 
important to understand the exact typ e of m a gnetic configurations wherein filaments 
form. One early classification scheme jTand. 119871 ) splits them into two categories 
based on the nature of the PIL above which the filament lies. The first category is 
one in which the filament forms above a PIL lying within a single bipolar unit of 
flux, and is classified as an "Internal Bipolar Region Filament" (see Figures 125a and 
126b ). In the second, the filament forms above a PIL which lies between two separate 
magnetic bipoles and is called an "External Bipolar Region Filament" (Figures l25b and 
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Fig. 25 Classification s cheme for solar filaments based on those of iTand (Il987|1 and 
Tandbcrg-Hansscn (1995) with two new categories (c and d) introduced by iMackav et al\ 
(200$). (a) Filaments which form above the internal PIL of a single bipole are classified as 
IBR. (b) Those forming on the external PIL between bipoles or between bipoles and unipolar 
regions of flux are classified as EBR. (c) Filaments that lie both above the internal PIL within 
a bipole and the external PIL outside the bipole are classified as I/EBR. (d) Filaments that 
form in diffuse bipolar distributions created through multiple flux emergences (such that the 
diffuse region can no longer be associated with any single bipole emergence) are classified as 
DB R, This category is ex pected to lie only at high latitudes. This figure is taken from Figure 
2 of lMackav et al\ d2008f> . 



125b). iTandberg-Hanssenl (|l995h descri bes th e se tw o categories as Type A and Type 
B, respectively. Early observations by iTand ( 19871 ) showed that, when filaments are 
classified into these two types, over 60% of filaments form external to bipolar regions. 

In a more recent study IMackav et all |200Sl) reconsidered where long, stable solar 
filaments form. Their study followed the history and evolution of the PILs underneath 
filaments using a wide range of data. The data included Ha synoptic maps, large- 
scale Ha images from the Ottawa River Solar Observatory (ORSO), Kitt Peak (KP) 
full disk and synoptic magnetograms, KP He 10830 syno ptic images, and results from 
magnetic flux transport simulations ( Yeates et all [2007). Through following the his- 
tory and evolution of the PILs above which filaments form, the authors determined the 
origin of flux regions at high latitudes from their initial bipolar source regions at active 
latitudes. As with previous studies the data and techniques employed were preferen- 
tially directed towards studying large, stable filaments, so they did not include smaller 
unstable filaments such as those found in the centers of activity complexes. 

To dis tinguish the d i fferen t bipole interactions that could lead to the formation of 
filaments. IMackav et al\ ( 2008J) introd uced two additional categories for filaments com- 
pared to those defined bv lTand (Il987l ): "Internal/External Bipolar Regions Filaments" 
(I/EBR) and "Diffuse Bipolar Region Filaments" (DBR). The I/EBR filaments are 
defined as filaments that lie above both the internal PIL of a bipole and the external 
PIL surrounding the bipole ( Figuresl25b a nd !26b ). and therefore could not be classified 
into the scheme proposed bv lTand (|1987l ). In contrast, the DBR filaments are located 
in essentially a bipolar distribution of flux, but where the polarities defining the bipole 
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Fig. 26 Examples of the four categories of filaments shown in Figure [25] (from M ackav et al\ . 
2008). In each of the panels (a)-(d), the bottom plot is an Ha image from the ORSO, while the 
top image shows the radial magnetic field derived from either (a)-(c) a full-disk magnctogram 
or (d) a synoptic magnetogram from Kitt Peak. Outlines of the Ha filaments are superimposed 
on each of the magnetograms. The dates of the observations are (a) 26*^ June 1979, (b) 6 th 
May 1979, (c) 27 th September 1979 and (d) 14 th July 1979. For panels (c) and (d) the areas 
enclosed by the boxes denote the corresponding area of (c) the magnetogram and (d) the Ha 
image. In panel (d) (top image) the low latitude activity complexes which will extend poleward 
over time and interact to produce diffuse regions of flux at high latitudes can be clearly seen. 



did not emerge together. The formation of the bipolar distribution was the result of 
many flux emergences, coalescences and cancellations such that the polarities on either 
side of the filament could not be attributed to a si ngle bipole eme rgence (Figures I25U 
and | 26£ . On comparing the clas sification scheme of lEngvoldl ( 19981 ) with those of lTand 
(|l987h and iMackav et ai\ i|200Sl ). ARFs could be of either Internal Bipolar, External 
Bipolar or Internal/External Bipole filament categories. IFs are always of External 
Bipolar type and QFs either Internal, External, Internal/External or Diffuse Bipolar 
Region type. Therefore in principle active region and quiescent filaments may form 
in very different magnetic environments. By using the f our categories (of IBR, EBR, 
I/EBR and DBR) and extending the work of lTand (|l987l ) to four distinct phases of the 
solar cycle (two before and two after cycle maximum) the authors were able to distin- 
guish more clearly the different types of bipole interactions leading to the formation of 
Intermediate and Quiescent filaments. 
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Of the 603 filaments studied by iMackav et al\ l|2008l ). 92% formed at locations 
requiring multiple bipole interactions (the breakdown comprised of 62% EBR, 17% 
DBR and 13% I/EBR). Only 7% formed within a single bipole. These results show 
that large-scale filaments, namely those of the IF and QF type preferentially form 
at sites of multiple bipole interactions. Very few of them occur within a single bipole . 
Furthermore, by considering four distinct phases of the solar cycle. IMackav et al . (2008) 
showed that only EBR filaments exhibit any form of solar cycle dependence, with th e 
other three types remaining essentially constant (see Figure 3 of IMackav et all I2008T ). 
The dependence showed that the number of EBR filaments varied in phase with the 
solar cycle with more at cycle maximum than minimum. Such a variation indicates 
that the formation of EBR filaments must be strongly related to the amount of flux 
emergence in the solar cycle. 



5.1.2 The Hemispheric Pattern of Solar Filament Channels and Filaments 

While the basic properties of solar filaments have long been known, fila ment channels 
and fi laments have been classified more recently in terms of their chirality |Martin et all 
hm4 ). As discussed in section 12.1.21 this chirality may take one of two forms: dextral 
or sinistral. Dextral/sinistral filament channels and filaments have an axial magnetic 
field that points to the right/left when the main axis of the filament is viewed from the 
positive polarity side of the PIL (see Figure [5J. The chirality of filament channels may 
be deduced from high resolution Ha images combined with magnetograms, whereas 
the chirality of filaments may be determined from that of the channel, direct mag- 
netic field measurements or from the relationship of filaments to their barbs (dextral 
filaments ~ right bearing barbs, sinistral filaments ~ left bearing barbs). In general 
due to the lack of high resolution Ha data and direct measurements of magnetic fields 
within prominences, filaments are most l y classified using th e ir rel ation ship to barbs 
Pevtsov ef"aH.l2003l;lYeates et all. 12007ft. lLopez Ariste et al\ (120061) and lMartin et al\ 



20081 ) used the chirality rules to resolve the 180-degree a mbiguity in photospheric 



vector-field measurements. In force-free field models (e.g., Ajilairiex aJid_rJemouhj, 
19981 ; IMackav. Longbottom fc Priestl . [l99 j|; lvan Ballegooiien et aZ.Tl2*0*0*oTlMackav*lind van BallegooiienI . 
2005) this chirality is directly related to the dominant sign of magnetic helicity that 
is contained within the filament and filament channel. A dextral filament will contain 
dominantly negative helicity, while a sinistral filament positive helicity. Hence, fila- 
ments and their channels may be regarded as indicators of sheared non-potential fields 
within the solar corona. Their transport across the solar surface is ther efore an indi- 
cation of the large-scale transport of magnetic helicity across the Sun l Yeates et al\ . 



2008lJ a key feature in explaining many eruptive phenomena. A surprising feature of 
the chirality of filaments is that it displays an unusual large-scale hemispheric pattern: 
dextral/sinistral filaments dominate in the northern /southern hemispheres respe ctively 
(jMartin et all Il994l ; Izirker et all I"l997l : IPevtsov et all I2OO3I : lYeates et adl2007l ). This 



pattern is unusual as it is exactly opposite to that expected from differential rotation 
acting on a North-South coronal arcade. Although dextral/sinistral filaments dominate 
in the northern/southern hemisphere, observations show that exceptions to this pat- 
tern do occur. Therefore any model which tries to explain the formation of filaments 
and their dominant axial magnetic fields must explain not only the origin of this hemi- 
spheric pattern but also why exceptions arise. The origin of this hemispheric pattern 
will be discussed in Section \5. 51 
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It is clear from the above discussion that solar filaments form or are found in a wide 
range of magnetic environments on the Sun, ranging from the rapidly evolving activity 
complexes to the slowly evolving poleward streams of flux that extend out of the active 
latitudes towards the poles. To explain the formation of these filaments, observational 
studies and a wide range of theoretical models have been produced. The review will 
now consider observational case studies of the formation of filaments (Section 15. 2|) . 
After discussing these, models of filaments formation will be discussed in Section [5.31 
The observations will then be used to clarify which models of filament formation are 
applicable to which filament formation locations ("Section 15. 6[) . 



5.2 Observations of Filament Channel and Filament Formation 

To understand the magnetic environment and interactions leading to the formation 
of filament channels and filaments, it is useful to discuss test cases. To date, very 
few examples of filament channel formation have ever been observed, so the exact 
formation mechanism remains debatable. Six recent publications present detailed case 
studies. Four cases show the formation of filament channels through surface effects that 
reconfigure pre-existing coronal fields, while in the latter two examples flux emergence 
of horizontal flux ropes is deduced by the authors to play a critical role. Thus from 
interpreting the observations there appear to be two opposing views on how filament 
channels and filaments form. In this discussion we will consider the key observational 
features and determine whether the two views may be reconciled. 



5.2.1 Evidence of the Reconfiguration of Pre-Existing Coronal Fields in the Formation 
of Filament Channels 



Observations reported by iGaizauskas et al. I (|l997h and IGaizauskas et al\ (|200lh show 



that surface motions acting on pre-existing coronal fields play a critical role in the 
formation of filament channels and filaments. In the first case, an Intermediate Filament 
(IF) forms over a short period of a few days, while in the second a Quiescent Filament 
(QF) forms over a period of months. In both cases the filaments form on PILs external 
to any single bipole and in the classification scheme of Section \5 . 1 . 1 1 would be classed as 
External Bipolar Regions Filaments. Although the two cases occur over very different 
time and length scales there are a number of important similarities. 

Both cases begin with the emergence of a significant amount of magnetic flux in the 
form of an activity complex. Importantly however, no filaments form during the process 
of flux emergence. In fact, for the large scale QF the filament forms approximately 27 
days after major flux emergence subsides. In both cases a necessary condition for the 
formation of the filament channels was flux convergence and cancellation at a PIL 
between separate bipolar regions. Such convergence and cancellati on of flux was also 
sho wn to b e impo rtant for filament formation in the papers by IGaizauskas! (2002) 



and iMarti n (1998). Finally and most importantly, in each case a significant amount 



of magnetic shear was seen to build up in the activity complexes as they emerged. 
The redistribution of this non-potential field through surface motions towards the PIL 
produces a preferred direction of the coronal field above the PIL and plays a critical 
role in the formation of the filament cha nnels, the necessary ingredient for filament 
formation ( Mackav and Gaizauskasl . I2003T I . 
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Fig. 27 He? images from Gaizausk as et al\ fiTlO V : (Figure 2) of the early stages in the for- 
mation of an Intermediate Filament between an old remnant region (bright plages A and 
B in bottom right images) and an emerging activity complex (inside oval). The Ha images 
correspond to 19 t,l -21 st July 1979 whereas the filament did not form until the 25 th July 1979. 



Figures |2"T1 and |2"51 illustrate the main stages in the formation of a filament channel 
and IF over a period of fiv e days between the 20 th -25 th July 1979 ( see Figures 2 and 4 
in lGaizauskas et all ll997T l. The formation of this southern- hemisphere IF involves the 
interaction of two distinct magnetic flux distributions, an old remnant region (M c Math 
16159) and a new emerging region (M c Math 16166). In the Ha image of Figure |2"T1 (top 
left) the bright North-South plage outlines the old remnant region; the oval denotes 
the location where new magnetic flux will emerge on the following day. A key feature 
of this image is that the chromosphere is free of any strong patterns of magnetic fields 
within the oval. New magnetic flux first emerges inside the rectangle within the oval 
on the 20 th July (bottom left) and then more strongly between the 20 th and 21 st (top 
right). Over this time there is very little change in the magnetic field of the old remnant 
region. Magnetic field observ ations show that the activity complex ( Gaizauskas et all , 
1 19831 ; iBenevolenskaval [2005) is made up of two or more sunspot pairs. Significantly, 
no filament forms near or around the activity complex during this period of rapid 
flux emergence. The key development in the formation of the filament channel occurs 
between the two right hand panels taken just three hours apart on the 21 s * July 
(at the location denoted by B in lower right panel). Over this period a band of co- 
aligned fibrils form at the tail end of the new activity complex, between it and the old 
remnant region. These co-aligned fibrils indicate a magnetic field at this location with a 
dominant horizontal component , i.e. that a filament channel has formed. According to 
a model bv lMackav et ali i 19971 ) this pattern of co-aligned fibrils can only be explained 
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Fig. 28 H« image (right) and magnetogram (left) from Gaizauskas et al. ( 199?) showing the 
final stages of the formation of an Intermediate Filament on the 24 th and 25 s * July 1979. For 
the magnetogram images white represents positive and black negative flux. The Intermediate 
Filament forms on the 25 th July after flux convergence and cancellation occurs at Fl. 



by the extended non-potential magnetic field of the activity complex in which the field 
contains a large amount of positive helicity (correct sign for the southern hemisphere). 



No filament forms as magnetic flux continues to emerge within the activity complex. 
The emergence ceases on the 23 rd July, thereafter the activity complex expands and flux 
begins to disperse into the filament channel. In Figure [28] the distribution of magnetic 
flux (left column) and corresponding Ha images (right column) can be seen for the 24 th 
and 25 th July. In the magnetogram images, white (black) represents positive (negative) 
flux. After major flux emergence ceases on the 23 rd the trailing positive polarity of the 
activity complex disperses or diffuses out. This dispersion causes a convergence of flux 
between the old and new regions. On 25 th July, five days after the complex started 
to emerge, cancellation of flux occurs at the point Fl (bottom right). The filament 
forms after this cancellation and passes through the location of flux cancellation. It 
was a stable structure which survived for a full solar rotation and can clearly be seen 
to lie on a PIL whi ch is external to any one bipolar region. Subsequent modeling by 
iMackav et all { 19971 ) showed that the resulting magnetic structure of the filament could 
only be explained by the interaction of the combined fields of both the old and new 
magnetic distributions. Both fields were highly non-potential, again with a significant 
amount of positive helicity which must have originated during the creation of the new 
activity complex. It is clear from the observations that reconfiguration of the previously 
emerged fields played a critical role in the formation of the filament channel. 
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Gaizauskas et al . (2001) described a similar process of filament channel and filament 
formation, but this time for a QF which is nearly IRq in length. The process of 
formation once again begins with the emergence of new flux, but this time in the form 
of two neighboring activity complexes in the northern hemisphere. Fibril alignment at 
the chromosphere once again shows that the extended fields of the activity complexes 
are highly non-potential with negative helicity. As with the previous case, the filament 
only forms after major flux emergence ceases and the activity complexes converge and 
partially cancel with one another. In contrast to the filament in the previous case, 
which took 5 days to form, for the large-scale case the QF takes nearly one full solar 
rotation (27 days) to appear. 

In both cases described above no stable filaments form during the periods of the 
highest rates of flux emergence, and the authors concluded that surface motions acting 
on pre-existing coronal fields play a critical role in the formation of stable filaments 
through the interaction of multiple bipoles. This result is consistent with the classifi- 
cation of filaments given in Section [5.1.11 where the majority of filaments are found to 
lie in magnetic configurations that involve more than one bipole. A key role of these 
surface motions is to redistribu te the helicity which is s e en to emerge in the early stages 
to fo r m the filament channel ([G aizausk as et all 1 19971 . l200ll : iMackav and Gaizauskas! . 
120031 ; I Mackav and van Ballegooiienl . I20051 . 120061 ). At a later stage a filament may ap- 
pear when mass is deposited into the channel. Convergence and cancellation of flux 
have also been shown to be important for filament channel formation according to 
the observa t ions r eported by iMartinl (|1998l ) who reviews several clear examples, and 
iGaizauskasI (|2002l ). who shows that early in the solar cycle a unipolar region of flux 
has to extend 180° around the Sun to interact and cancel with an opposite polarity 
region before a filament can form on that PIL. For these cases the redistribution of 
flux, after emergence, is inferred to b e a key pro c ess in the formation of the filament 
channel. In the more recent paper bv IGaizauskasI d2008T ). three different arrangements 
of interacting nests are considered for the formation of filament channels and filaments 
as flux disperses across the Sun. In each case channels and quiescent filaments form on 
the boundaries of these activity nests. 

ISchmieder et al] (|2004h studied the formation of a filament in the complex cen- 
ter of a decaying active region formed out of smaller individual components, using 
multi-wavelength observations obtained during a "Joint Observing Programme" be- 
tween ground-based instruments in the Canary Islands (the SVST and the MSDP on 
the VTT) and the TRACE satellite. They followed the evolution of three individual 
filament segments denoted Fl, F2 and F3 over several days, and found that Fl and 
F2 gently merged into a single structure, as observed by a gradual filling in Ha of the 
gap between them. This merging was associated with mild EUV brightenings and with 
small Ha Doppler shifts at the merging point. While EUV brightenings are a good 
indicator of magnetic reconnection, the flows revealed that the merging first took place 
by dynamic exchanges between the two progenitors, until they formed a single long 
stable filament. Two days later segments F2 and F3 c ame into contact and produced a 
confined flare, as evidenced by EUV post-flare loops (foeng et a/.l. l2002ll . To determine 
the directions of the axial fields in the three filament segments, Schmieder et al\ |2004) 
used the chirality rules for chromospheric fibrils and magnetic field polarity, the skew 
of the overlying coronal arcades, and the sense of twist in neighboring sunspots. It was 
then confirmed that when two filaments interact, magnetic reconnection takes place 
and leads to a merging when their chiralit ies are of the same sign, but leads to a flare 
when the chiralities are opposite (also see iMalherbel. Il989l ; IMartinl . Il998l ; iRustl l200ll ; 



53 



van Ballegooiien1 . l2004l ; lDe"Vore et <afl l2005h . It was also inferred that magnetic helicity 
must slowly accumulate prior to filament merging, as seen by the rotation of a small 
twisted sunspot close to the merging point. Finally, it was suggested that magnetic 
reconnection first accelerates plasma between both progenitor filaments, and that it 
may later result in a change of topology which can sustain stable plasma all along the 
new filament. 

I 1 J 

More recent observations by Wang and Muglach (2007) have supported the work 



of iGaizauskas et al. | (|l997l , bOQlMWang and Muglacbl (|2007l ) describe examples of the 



formation of filament channels and filaments through comparing BBSO Ha images with 
MDI normal component magnetograms. The authors describe how fibrils which are ini- 
tially normal to the PIL rotate to lie parallel to the PIL over a period of I to 2 days and 
in doing so form a filament channel. Through studying the evolution of the magnetic 
fields the authors deduce that flux cancellation as a result of supergranular convection 
plays a key role in the formation of the filament channels. They argue that this can- 
cellation process between opposite polarity elements removes the normal component 
of the field but leaves the component parallel to the PIL w hich builds up gradually 
to fo rm the axial field of the filament channel. In contrast to IGaizauskas et al\ (|l997l . 
l200ll l they do not observe any significant helicity resulting from the emergence of the 
active regions. For the two clearest examples of filame nt formation, 15th January 2002 
and October 7th 2002 (see IWang and Muglachl l|2007l ') Fi gures 2 and 5) the filaments 
form between bipolar regions of flux and background fields and therefore would be 
classified as IF and Exterior Bipolar Region Filaments. For the third case which is not 
so clear, the filament partially lies on both the internal and also the external portions 

of the PIL so would be an I/E BR filament or an AR fi lament. 

On comparing the results o flGaizauskas et al\ ( 19971 ) and lWang and Muglachl ( 20071 ') , 



while there are many similarities, there are also some differences in time scale. The 
cleare st is the time difference required to form the filament channel. For lGaizauskas et al\ 
(|l997T ) the formation of the filament channel occurs over a 3 hour period and is at- 
tributed to the extended non-potential field of the activity complex containing a large 
amount of helicity. Cancellation of flux could not produce such a strongl y sheared field 
over such a short period of time. In contrast. I Wang and Muglachl (2007) do not report 
any strong patterns of fibrils associated with helicity emerging in the active regions but 
rather form the filament channel over a period of 1 to 2 days in a much slower process 
of cancellation. Therefore there appear to be two complementary methods of forming 
a filament channel over different time scales. 

5.2.2 Evidence of Emerging Horizontal Flux Tubes in Filament Formation 

It is clear from the above observations that surface effects play a critical role in 
formin g the studied IFs and QFs (which are long stable structures). iLites and Low! 



(1997) describe a differe nt process for forming short, unstable active-region filaments. 
In lLites and Low! l|l997h the emergence of a 5-spot is traced through vector magnetic 
field measurements using Advanced Stokes Polarimetry. Magnetic field vectors along 
part of the PIL within t he emerging £- s pot sh ow a co ncave up or dipped magnetic 
structure (see Figure 1 of ILites and Low! 1 19971 ); also see ILites] (|2005l )). A small active 



region filament forms at this locatio n . The filament was however unstable with a life- 
time of only 2 days. ILites and Low! |l997) suggest photospheric material is dragged 
up into the corona thro ugh the levitation process, as a horizontal flux rope emerges 
( Rus t and Kumarl <| 1994b , also see section [3]). 
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A more recent example of the effect of evolving mag netic fields on the stru cture 
and stability of an ac tive region filament is described by lOkamoto et al\ ((20081 ) and 
lOkamoto et al.l ( 20091 ). In two papers, the authors present observations of a time series 
of vector magnetic fields taken by SOT underneath a pre-existing filament. The vector 
magnetic field measurements show a PIL with dominant horizontal field along it. This 
horizontal field probably represents that of the filament channel of the pre-existing 
filament. Over a period of 1.5 days the horizontal field vector changes from normal to 
inverse polarity and a dominant blue shift is observed. During this period the filament 
alters its appearance from a single structure, to a fragmented one and back again. 
Before returning to a single structure, brightenings are observed along the filament 
fragments in the Ca II H line. 

From the observations the authors deduce two possible scenarios. In the first sce- 
nario they interprete the observations in terms of an emerging horizontal flux rope 
which fully emerges into the corona and occupies the position of the pre-existing fila- 
ment. With this scenario the mass of the prominence originates from below the pho- 
tosphere. The second scenario interprets the brightenings in Ca II H as evidence for 
reconnection between the pre-existing filament and a new flux rope that emerges free 
of mass. The reconnection then produces a single structure along the PIL. A difficulty 
with both scenarios is that no simulations of magnetic flux emergence have been able 
to emerge a horizontal flux rope through the photosphere. 

In contrast too that put forward by the authors, a third possibility also exists. 
This fits the theoretical models discussed in Section T3. 21 As the top part of a flux rope 
emerges, a likely outcome is the emergence of sheared arcades. A coronal flux rope 
may then be formed out of these arcades through the process of reconnection. This 
reconnection may lift cool material into the corona, as has been discussed in Section [3] 
If the axial component of the emerging arcade lies in the same direction as that of 
the pre-existing filament channel, the new and old flux systems may join to produce a 
single structure (see Section[5]4|. To consider which, if any of these three scenarios are 
correct, new high resolution magnetic field observations at multiple levels in the solar 
atmosphere (e.g. photosphere, chromosphere and corona) are required. 



5.2.3 Summary of Observations 



The observations described in Section 15.2.1 1 and 15.2.21 provide evidence for filament 
formation arising from surface motions that reconfigure already existing coronal fields 
or, emerging flux tubes. So can the two methods be reconciled? The important dis- 
tinction between these cases is the type and location of filaments formed in each case. 
For the first four cases surface motions play an important role in forming long stable 
Quiescent or Intermediate filaments which are External Bipolar Region Filaments, the 
dominant type of large-scale filament found at all latitudes on the Sun. In contrast, 
flux tubes emerging in a <5-spot forms an Active Region or Internal Bipolar Region 
Filament which is unstable, lasting merely two days. 

While it is difficult to draw general conclusions from just six specific observations, 
they indicate that two different mechanisms might form filaments in different magnetic 
environments on the Sun. Thus large stable filaments of the IF and QF type (External 
or Diffuse Bipolar Region) may require surface motions to gradually reconfigure pre- 
existing coronal fields, while small, short-lived ARFs (Internal Bipolar Region) may 
form due to flux emergence. To determine whether different mechanisms do produce 
different types of filaments at different locations on the Sun, the formation of filaments 
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Table 3 Surface Models of Filament Formation 



Single Bipolc 

van Ballegooiien and Martens (1989) 1 ' 3 ' 4 ' 1 ' 
DeVore and Antiochos (2000) 1 ^~ 



Multiple Bipoles 



Kuperus (1996) 1 ' 3 - 4 

Kuijpers q9971 3 ' 4 - 8 ' 10 

Mackav et al. (19981 3 ' 4 - 6 ' 8 - 10 

Galsgaard and Longbottom (^99) 3,4 

van Ballegooiien et al. (20001 1 ' 4 ' 10 

Martens and Zwaan (2001 i 3 - 4 ' 10 

Lionello et al. (2002) 8 - lu 

DeVore et al. (2005) 1 - 3 ' 4 

Mackav and van Ballegooiien (2005) 1,4,8, 

Welsch et al. (2005) 3 ^' 8 - 10 

Litvinenko and Wheatland (2005) 3 ' 4 ' 8 ' 10 

Yeates et al. ('2008a') 1 ' 4 - 8 ' 10 



Table 4 Sub-Surface Models of Filament Formation 



Single Bipolc 



Multiple Bipoles 



Low (19941 7 

Rust and Kumar (1994) 7 ' 9 
Gibson et al. (2004V ''' M 
JjOW_and Hundhausen M995) 7 ' 9 
Fan and Gibson (2004l 7 ^ 
Fan and Gibson (2006) 7 ' 9 
Gibson and Fan (2006) 7 ' 9 
Magara (2006) 7 - 9 
Magara et al. C2008) 7 - 9 
(Fan, 200916. 9 



van Ballegooiien and Martens (1990) 2,3 ' 4 ' 7 
Priest et al. (WdG) ' 2 ^'^ 
Oliver et al. (1999) 2 ' 3 ' 4 ' 6 



over a wide range of latitudes needs to be considered in detail. Observational programs 
required to do this will be discussed in Section 15.6.31 



5.3 Theoretical Models of Filament Formation 

Over the years many models have been constructed, each employing a variety of mech- 
anisms in order to describe the formation of filaments. These models vary from descrip- 
tive papers to full numerical MHD simulations and consider two main problems. First, 
how to obtain the correct dipped magnetic field configuration with dominant axial mag- 
netic field that follows the hemispheric pattern, and secondly, the origi n of the dense 



plasm a. While the second question relates more to thermodynamics (jKarpen et al\ 
(2001), section[3]), this section which relates filaments to their underlying magnetic po- 
larities is relevant to the first group of models. It is widely accepted that magnetic flux 
ropes are a suitable configuration to represent solar filaments; the main area of debate 
is how exactly these flux ropes may form. Therefore, the various models in that group 
may be broadly split into two distinct sub-groups: those employing surface effects to 
reconfigure coronal fields (Table [3j) and those employing subsurface effects (Table [4]). 
This split naturally arises from the discussion of the observations in Section 15.21 At 
the present time only those employing surface mechanisms can be directly compared to 
observations. In these tables the surface/subsurface models have also been subdivided 
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Table 5 Mechanisms of Filament Formation 



Surface Mechanisms 

(1) Differential Rotation (shear flows) 

(3) Converging Flows 

(4) Magnetic Reconnection (atmosphere) 
(6) Flux Emergence (bipoles) 

(8) Magnetic Helicity 

(fO) Flux Cancellation/Diffusion 



Subsurface Mechanisms 

(2) Subsurface Motions 

(5) Magnetic Reconnection (subsurface) 
(7) Flux Emergence (U-loops) 
(9) Magnetic Helicity 



into those acting in single or multiple bipolar configurations in account of the observa- 
tions discussed in Section [5.1.11 The list contains models which consider the physical 
processes and mechanisms that may produce the 3D magnetic structure of filaments. It 
should only be regarded as representative and not exhaustive. Where the same authors 
publish multiple papers based on the same mechanism, generally only the first paper 
outlining the process is listed. Due to this, readers are recommend to search for other 
such papers in the literature. For each of the entries in the table the numbers attached 
correspond to the various mechanisms that the models employ, as listed in Table [5] 
From Table [3] it is clear that surface models rely on a variety of mechanisms combined 
together, while subsurface models generally rely on the emergence of twisted flux ropes 
where the filament forms in the dips of the flux rope or U-loop. 

Surface mechanisms include: differential rotation; shear flows along a PIL (differ- 
ential rotation is just a we ak shear flow); and co nverging flows onto a PIL. Recent 
helioseismic observations bv lHindman et all (|2006h show that underneath a well devel- 
oped filament strong shear flows may be observed. However this was after the filament 
had formed and not during the formation process. Photos pheric converging or shearing 
flows may be detected by local correlation tracking (cf. iMaearal . 1 19991 : iRondi et all . 
2007; Roudier et all 12008). For some surface models diffusion of flux towards a PIL 
with subsequent cancellation plays the role of the converging flow. To produce an axial 
magnetic field consistent with observations, these flows generally have to occur in a 
specific order. In contrast to these surface motions, subsurface shear flows have also 
been employed. In both sets of models magnetic reconnection is generally required to 
reconfigure the fields; the reconnection may occur either above or below the surface. 

Another feature common to both sets of models is flux emergence, but it is used in 
very different ways. For surface models, magnetic bipoles which emerge either untwisted 
or twisted are advected across the solar surface and reconfigured with other pre-existing 
coronal fields as discussed in the observations of Section A key element in recent 

papers describing filament fo rmation is that these bipoles are non-potential and includ e 
an initial magnetic helicity (Mack av and van BallegooiienL 120051 : lYeates et all l2008al ) . 
In contrast, flux emergence for subsurface models is presumed to occur in the form of 
twisted U-loops (Section 15.2.21) . 

Whilst it is impractical to describe each of the models listed in Tables [3] and [4] in 
detail, key elements may be considered from a few selected cases. The cases chosen 
are picked solely for illustr ative purpos es. T he key feature of any sub -surface model is 
described in the papers bv lLowl l|l994 ) and iRust and Kumarl (l994 ). For these mod- 
els a filament is formed by a horizontal twisted magnetic flux tube in the convective 
layer of the Sun. Due to magnetic buoyancy the tube rises through the convective 
layer and emerges into the photosphere and corona, dragging cool dense material with 
it, to produce the mass of the prominence. Not every part of the tube has to rise at 
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the same ti me; subsequent ris i ng; pa rts could explain canceling magnetic features (see 
Figure 13 in lRust and Kuman. 19941). Such a fe ature has been consid ered in the numer- 
ical si mulations of lGibso'ri*lmdF=mf ( 20061 ) and iMagara et al\ ( 20081 ) . In iMagara et al\ 
(2008) the flux rope is forced to rise into the solar atmosphere by imposed velocity 
fields. In other flux emergence simulations where the authors use only buoyancy and 
magnetic buoyancy instabilit ies, it is found that the axis of the flux rope does not 
rise through the photosph e re jMoreno-Insertis , 20041 ; lArchontis et aUl2004l ; [Archontisl . 
120081 ; iMurrav et aH |2006| ; IGalsgaard et all 120071 ) Although the axis and U-loops of 
the emerging tube do not rise to coronal heights, the process of flux emergence may still 
produce a coronal flux rope with dips. A flux rope may form throug h the reconnection 
of emerged sheared field lines that lie above the emerging tub es axis ^Manchester et all . 

120041 ; iMagaral . |2006| ; lArchontis and TorokL |2008|; iFarl [2009I). 

In contrast, one of the first surface models, bv lvan Ballegooiien and Martens! ( 19891 ). 
considers shearing motions acting on a coronal arcade in a bipolar configuration. The 
footpoints of the arcade are sheared in such a way that their separation increases and an 
axial field component is produced along the PIL (see Figures [29k and !29b ). In principle 
this shear could be supplied by solar differential rotation or by other shear flows on 
the Sun. Next convergence, or diffusion of the flux towards the PIL, brings the foot 
points together where they may reconnect to produce, a long axial field line along the 
PIL and also a small loop which submerges through the surface (Figures I29fc and 1291 1). 
Subsequent repetition of this process creates dipped magnetic field lines consistent with 
the topology required for filaments (Figures l2"9"fe andl2"9"f ) 



In an extension to the lvan Ballegooiien and Mar tens ( 1989]) model. lMartens and Zwaanl 



(20011) put forward a "head-to-tail" linkage model for the formation of filaments in a 
multiple bipolar configuration. A key element of their model was once again flux con- 
vergence and cancellation, which acts a s the driver for reconnect i on be tween initially 
unco nnected magnetic sources ( also see iKuperusl 1 19961 : iKuiipersl Il997l )). In this sce- 
nario jMartras_^ncT^waarj ( 200 ll ) describe how the filament channel and filament may 
be produced by the interaction of two bipoles, one older and more diffuse lying at 
a slightly higher latitude. As long as these sources satisfy Hale's Polarity Law and 
Joy's Law, then convergence and cancellation as a consequence of differential rota- 
tion, of the following polarity of the lower latitude bipole and lead polarity of the 
high latitude bipole, could result i n a st rongly sheared field line along the PIL. As 
with Ivan" Ballegooii en and Martens] ( 19891 ) successive repetitions of this process would 
build up helical field structures. In addition, if multiple bipoles are involved over a 
range of latitudes large filaments extending over a full solar radius could be produced. 
While Martens and Zwaan conside r ed the interaction of multiple bipol es in a con- 
ceptual model, iMackav et al\ (|2000l ) ; iMackav and van Ballegooiienl (|200ll ) carried out 
numerical simulations of a similar process. 

IGalsgaard and LongbottomI ( 19991 ) consider a similar scenario with 3D numerical 
MHD simulations without the processes of flux cancellation and submergence. The 
interaction of two bipolar regions of flux, one which is older and more spread out 
than the other, is considere d. The numerical s imula tions are based on the magnetic 
configurations discussed by iGaizauskas et al. I jl997l) . In the model the two bipoles 
are initially connected to one another and the authors demonstrate how a current 
sheet may form in the coronal volume between the two bipoles as a result of flux 
convergence but no cancellation or submergence. Subsequent reconnections resulting 
from the convergence, then lift cool dense plasma over a number of pressure scale 
heights where it is able to form a region of high density plasma overlying the PIL. As 
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(d) (•) tfl 

Fig. 29 Example of the formation of a filament's axial magn etic field through shearing mo - 
tions, convergence, and reconnection as put forward by van Ballcgooiicn and Martens (1989). 



part of this reconnection process dipped and sometimes helical field lines are formed 
with an upward tension force. Although the authors carried out the simulation with 
a simple isothermal atmosphere it remains one of the few full MHD simulations to 
consider both the origin of mass and of the shear in the magnetic field. 

An alternative method o f forming a similar magnetic structure was proposed by 
DeVo re and Antiochosl (|2000l ) using a single bipolar configuration (also see lAntiochos et all 



19941 ). In this model, a bipolar magnetic field distribution is subjected to a strong 



shearing motion parallel to the PIL, however, no converging flow is applied. Once the 
footpoints of the field lines are sheared a distance comparable to the bipole width, an 
untwisted dipped magnetic configuration forms. The authors show that through fur- 
ther shearing of the dipped field lines the initially untwis t ed fie ld may form a helical 
structure similar to that of Ivan Ballegooiien and Martens! (ll989Tl through a tw o stage 



reco nnection process. Therefor e, in contrast to van Ballcgooi ien and Martens! (1989) 
and iMartens and Zwaanl ( 200 ll ) , iDeVore and Antio^hosTfeoOOl ) do not rely on conver- 
gence and cancellation of flux to produce the helical field. 

From the discussion above it is clear that a wide range of theoretical models exist to 
explain the 3D magnetic structure of solar filaments. These models employ a variety of 
mechanisms. As will be discussed in Section l531 at the present time none of these models 
may be ruled out. However, by combining the observations discussed in Sections 15.11 
and 15.21 it may be argued that some models are more relevant than others for the 
formation of large stable filaments (Quiescent and Intermediate) compared to Active 
Region filaments. A full discussion along with the presented hypothesis will be carried 
out in Section T5, 61 



5.4 Models of Filament Merging 

Numerical MHD simulati ons of the formation and interaction between filament seg- 
ments were conducted bv lDeVore et al\ ( 20051 ) . The footpoint motions of the sheared 
arcade model (section 12.2.21) were applied to two adjacent, initially current-free, mag- 
netic dipoles. As shown in Fig. 1301 four possible combinations of chiralities (identical 
or opposite) and axial magnetic fields (aligned or opposed) between the participat- 
ing filaments were considered. These four simulations exhibited substantially different 
degrees of linkage. 
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Fig. 30 Schematic diagram of the four possible configurations of interacting prominences: 
identical (top) or opposite (bottom) chiralities, and aligned (left) or opposed (right) axial 
magnetic fields. Black lines are the polarity inversion lines of the vertical field, whose direction 
is denoted by + (upward) and — (downward); shaded gray rectangles are the prominences, 
whose chiralities are indicated by d (dextral) and s (sinistral) and whose axial field directions 
are shown by the arrows; and filled black circles mark the prominence interaction regions. 



When a single polarity inversion line is shared by the two bipoles, then identical 
chiralities necessarily imply aligned axial fields (Fig. 1301 left). In this case, magnetic 
reconnection between field lines from both bipoles links the two initial prominence seg- 
ments. Both acoustic and Alfven waves propagate along these newly reconnected field 
lines, and should be capable of driving existing plasma condensations from one pro- 
genitor to another. As the shear increases, the volume between the original segments 
becomes filled with reconnected field lines, so that they gradually merge into a sin- 
gle filament. This multi-step merging mechanism couples photospheric shear, coronal 
reconnection, and relaxation to form a longer st ructure contain i ng nu merous dipped 
field lines where plasma can collect most easily ijAulanier et al\ . 120 06). Furthermore, 
the case of identical chiralities and aligned axial fields success f ully reproduces the ob- 
servations of filament merging (e.g., iMartin et all. 1 19941 ; iRustj l200ll ; | Schmieder et all 
2004; Ivan Ballegooiienl .[2004). A second model configuration (Fig. 1301 upper right) also 
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contains a shared PIL, but opposite helicities and axial fields between the two segments. 
No merging ensues because little reconn ection occurs betwe e n the filament segments , 
as is consistent with observ ations (e.g., iMartin et all 1 19941 ; iRustl . 1200 ll ; iDeng et~al\ , 
120021 ; ISchmieder et aUl2oblb . 

When the initial topology instead is quadrupolar, so that the system contains two 
orthogonal PILs (Fig. 1301 lower panels), then the converse relation holds between chiral- 
ity and axial-field alignment. Reconnections that form linking field lines now occur only 
between filament segments with opposite chiralities (Fig. 1301 lower right), while recon- 
nection between same-chirality segments only results in footpoint exchanges (Fig. 1301 
lower left). The results for quadrupolar topologies have not yet been verified with solar 
data because multipolar filament interactions are rarely observed. However, these key 
predictions present an important objective for future observational campaigns. 



5.5 Origin of the Hemispheric Pattern of Filaments 

Any model which tries to explain the origin of a filament's magnetic field must also ex- 
plain why this magnetic field exhibits a hemispheric pattern. The first attempt th rough 
detailed numerical simulations was carried out by Ivan Ballegooiien et al\ l|l998b who 
considered whether the surface flux transport effects of differential rotation, meridional 
flows and supergranular diffusion could in fact create the observed axial fields in fila- 
ments. By using observed magnetic flux distributions and initial coronal fields which 
were potential, they simulated the evolution of both the photospheric and coronal 
fields. They found that the above-mentioned surface effects when acting on potential 
fields create approximately equal numbers of dextral and sinistral channels in each 
hemisphere, in contradiction with observations. However, these authors did not take 
into account the force balance of the coronal plasma, and only considered statistical 

relationships between filament c hirality and latitude. 

In the more recent study of iMackav and van Ballegooiier] l|2005h the authors re- 



consider the origin of the hemispheric pattern through combined flux transport and 
magneto-frictional relaxation simulations. In these simulations as the flux transport 
effects shear up the photospheric field, the coronal field responds to these motions 
by relaxing to a non-linear force-free field equilibrium. In the simulations the authors 
consider an idealized setup of two initially isolated interacting bipoles. They perform 
a parameter study to determine the type of chirality formed along the PIL between 
the two bipoles, as the bipoles initia l tilt angle (Joy's law) and h elicity are varied. 
Therefore in contrast to the study of van Ballegooii en et al\ l|l998b they did not use 
initial potential fields. Through the simulations the authors demonstrate that surface 
diffusion can play a key role in canceling flux between the bipoles and building up 
an axial field through transporting helicity and sheared fields from the inner parts of 
the bipoles to the outer edges. The authors demonstrate that the hemispheric pattern 
of filaments may be explained through the observatio nal properties of newly em erging 



bipoles such as, (1) their dominant tilt angles (-10:30°. IWang and Sheeley 1 19891)). and 



(2) t he dominant helicity that they emerge with in each hemisphere (|Pevtsov et all 
1995). In addition to this, a key feature of the simulations was that for the first time 
the occurrence of exceptions to the hemispheric pattern as a result of large positive 
bipole tilt angles and minority helicity could be quant ified. 

The results of lMackav and van" Ballegooiien! J2005h have been tested throu gh a di- 



rect comparison between theory and observations bv lYeates eFafl ( |2007l[2008e^ . In this 
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Fig. 3 1 Example of the comparison of theory and observations performed by lYeates et al\ 
(2008a). (a) Example of the magnetic field distribution in the global simulation after 109 days 
of evolution, showing highly twisted flux ropes, weakly sheared arcades and near potential 
open fields. On the grey-scale image white/black represents positive/negative flux, (b) Close 
up view of a dextral flux rope lying above a PIL within the simulation, (c) BBSO Ha image 
of the dextral filament observed at this location. 



comparison lYeates et al. I l|2007l , l2008al l develop a new technique to model the long term, 
global evolution of the Suns magnetic field based on act ual magnetogram observations. 
To carry out the comparison lYeates et al ] (|2007l . l2008ah first consider Ha observations 
from BBS O over a 6 month period and determine the location and chirality of 109 
filaments ( Yeates et all 120071 ) relative to the underlying magnetic flux. It should be 
noted that all of the filaments in the study lie below 65° latitude and no polar crown 
filaments were included since their chirality could not be determined. To determine 
the chirality for such filaments direct magnetic field measurements would be required. 
In the second stage they use combined magnetic flux transport and magneto-frictional 
relaxation simulations, where the simulations are based on actual photospheric mag- 
netic distributions found on the sun. They show that by including the process of flux 
emergence in the flux transport simulations they can reproduce to a high accuracy KP 
synoptic magnetograms from one Carrington Rotation to the next. Over the 6 month 
simulation to maintain the accuracy, 119 bipoles were emerged with their properties 
determined from observations. When the bipoles were emerged they were initially iso- 
lated from the surrounding fields and could contain either positive or negative helicity. 
Their subsequent interaction with neighboring field regions and transport across the 
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solar surface could then be followed. A key difference of these simulations from previous 
studies, was that they were run for the whole six month period without ever resetting 
the surface field back to that found in observations or the coronal field to potential. 
Therefore the simulations were able to consider long term helicity transport across the 
solar surface from low to high latitudes. 

Using this technique I Yeates el al. (2008a) carried out a direct one-to-one compari- 
son of the chirality produced by the model with the observed chirality of the filaments 
at the exact location that the filaments were observed. An example of this can be seen 
in Figure 1311 In Figure 131b the global field distribution can be seen after 108 days of 
evolution. The photospheric field distribution is given by the grey-scale image (white 
is positive flux, black is negative flux), while the lines denote the field lines of the 
non-linear force-free coronal field. 

The coronal field can be seen to be made up of highly twisted flux ropes, slightly 
sheared coronal arcades and near potential open field lines. A zoomed in area of the 
central portion of this image can be seen in Figure [31b . In this figure a simulated flux 
rope structure can be seen where the axial field is of dextral chirality. For comparison 
the Ha filament that formed at this location can be seen in Figure 131b . Through 
studying the barbs and applying a statistical technique, the filament was determined 
to be of dextral chirality so the chirality formed in the simulation matches that of the 
filament. 

Through varyin g the sign and amount of helicity emerging within the bipoles, 
lYeates el al\ 1 2008al ) (see their Figure 5b) show that by emerging dominantly negative 
helicity in the northern hemisphere and positive in the southern, a 96% agreement 
can be found between the observations and simulations where the agreement is equally 
good for minority chirality filaments as well as for dominant chirality filaments. A 
key feature of the simulations is that a better agreement between the observations and 
simulations is found the longer the simulations are run. This indicates that the Sun has 
a long term memory of the transport of helicity from low to high latitudes. The reason 
for this high agreement is described in the paper of I Yeates and Mackavl 1 20091 ) where 
seven different mechanisms are involved in producing the observed chiralities. The 
results demonstrate that the combined effects of differential rotation, meridional flow, 
supergranular diffusion and emerging twisted active regions are sufficient to produce 
the observed hemispheric pattern of filaments. 



5.6 Discussion of Filament Classification Schemes, Observations of Formation and 
Theoretical Models 

So far we have discussed a number of properties of solar filaments. These properties 
have ranged from observations and locations of filament formation, to the wide variety 
of theoretical models used to explain them. In this sub-section we turn our attention to 
tying all of these observations together, by forming a unifying hypothesis, to quantify 
where and at what locations the mechanisms and models discussed in Section [5.31 are 
appropriate. We hope that this hypothesis will stimulate new observational studies 
to test it. We will first discuss the formation of large-stable filaments, namely those 
of the Intermediate and Quiescent filament type. After this we will consider Active 
Region filaments. A distinction is made as present observational evidence suggests 
that different mechanisms may apply. 
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5. 6. 1 Intermediate and Quiescent Filament Formation 

From the observations of large-scale stable filaments discussed in Section r5.1.1l it can be 
seen that the majority of IF and QF (92%) form in magnetic configurations involving 
multiple bipole interactions. Very few (7%) form in single bipolar configurations along 
the internal PIL of the bipole. From this it is clear that it is rare to find either IF or 
QF in such single bipole configurations. Therefore, while none of the models listed in 
Table0or[3]can be ruled out, it is clear that those involving multiple bipole interactions 
are the most appropriate. 

The question now turns to whether the IF and QF are formed d ue to surface mo- 
tions actin g on pre-existing coronal fie lds as lGaizauskas et al\ (jl997T l. IGaizauskas et all 
(|200lh and lWang and MuglachI (|2007|^ argue (Sect ion 15.2.1 [I or whether they a re due to 



sub-surface processes as argued bv lLites and Low! ( 1997h andlOkamoto et all J2009l 



The answer to this may be alluded from the observations of Mac kav et al.\ (2008) 
who showed (see their Figure 4) that only External Bipolar Region Filaments (EBR), 
those forming between two neighboring bipoles show any form of solar cycle depen- 
dence. Their numbers were found to increase in phase with the solar cycle, peaking 
around solar maximum and with more present in the declining phase of the cycle than 
in the rising phase. In contrast, Internal Bipolar Region Filaments (IBR) those form- 
ing in a single bipole, clearly showed no solar cycle dependence, with their number 
remaining at a low constant value throughout the cycle. If flux emergence in the form 
of flux ropes as used by subsurface models is key to the formation of IF and QF, then 
one would expect that IBR Filaments should show a strong solar cycle dependence. In 
addition they should be the dominant type. As more bipoles emerge on the Sun there 
is a much greater chance of having emerging flux ropes, hence more Internal Bipolar 
Region Filaments. As they are not the dominant type and exhibit no solar cycle de- 
pendence, another process other than flux rope emergence must be acting to form IBR 
filaments. An alternative explanation for the formation of IBR filaments, which from 
the observations discussed in this review generally occurs outside activity complexes, 
is that they may result from strong surface shearing motions or helicity injection and 
flux cancellation. Evidence for such shearing motions, at least i n the later development 
of a filament, has recently been published by Hindm an et al This fits in with 



the single bipole model of iDeVore and Antiochos! (|2000h where shear flows play a key 
role in the formation of the filament. In contrast, the alternate process of injectio n of 
helicity followed by cancellation has been described by I Wang and MuglachI (|2007h . 

If flux rope emergence is unlikely to form IBR filaments, then it is extremely unlikely 
that such a process could apply to External Bipolar Region Filaments (EBR), the 
dominant type which tend to form between distinct bipoles after flux emergence has 
ended (see Section 15.2.11) . Again, another process must be acting, one that is closely 

related to the amount of magnetic flux on the Sun. 

The obvious c hoice dedu ced from the observations of iGaizauskas et al. 

I (| 19971 ). 



IGaizauskas et al\ (|200lh and IWang and MuglachI (|2007r ) is the convergence between 
individual bipoles resulting in flux cancellation and reconnection, since the rates of 
convergence increase due to the increased rate of flux emergence during periods of high 
activity. During these periods of high activity the most widespread source of conver- 
gent flows would be the natural expansion of all bipolar pairs as soon as they emerge. 
If these processes are key in explaining the formation of filaments this would lead to 
an increase in the number of filaments forming between individual magnetic bipoles as 
activity increases. In addition, one expects more EBR Filaments after cycle maximum 
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than after cycle minimum; after maximum there is still a significant amount of mag- 
netic flux on the Sun so that convergence and cancellation or magnetic reconnection 
can still take place. 

Therefore, it may be argued that convergence leading to subsequent cancellation 
and reconnection (i.e., items 3, 4 and 10 in Table[S]) are the mechanisms that result in 
the formation of the majority of large stable filaments found on the Sun, and that flux 
rope emergence does not play a major role for these filaments. The models in Table [3] 
which include these mechanisms appear to be the most appropriate. At the present 
time no further distinction can be made between these models. While this argument 
is put forward for EBR Filaments, it also applies to Diffuse Bipolar Region Filaments 
(DBR) because the flux patterns in which they form are a natural consequence of flux 
convergence and cancellation occurring over long periods of time. 

Finally, for Internal/External Bipolar Region Filaments again convergence and can- 
cellation/reconnection may apply but this time strong shear flows or helicity injection 
may again be applicable. While the above discussion applies to Intermediate and Qui- 
escent filaments, in the next sub-section we consider Active Region filaments. 



5.6.2 Active Region Filaments 

We now turn our attention to Active Reg i on fil amen ts to consider their f ormation 



mechanism. Observations bv lLites and Low 1 19971 ) and lOkamoto et al.l ( 20091 ) suggest 



that small-scale unstable active region filaments may be formed as the result of flux 
rope emergence dragging cool dense photospheric plasma into the corona. While this 
remains a possibility, most numerical simulations of emerging flux ropes fail to lift the 
axis and cool material of the original flux tube into the corona. Therefore is remains 
unclear whether such a process may occur. In contrast, many authors have shown 
that during the process of flux emergence, a fter the top of the flux ro p e has emerged , 
magnetic reconnection or helicity inj ection ^Manchester et aZl |2004| ; iMagaral . 120061 : 
lArchontis and Torokl . 120081 ; iFanl . l2009f) may reconfigure the emerged coronal arcade to 
produce a secondary coronal flux rope. During the formation of the secondary flux rope 
the reconnection may then lift cool dense material to coronal heights. Therefore while 
emerging flux appears to be important for the formation of active region filaments a 
key element may still be atmospheric reconnection of pre-emerged fields. To resolve 
this issue many new observational studies are required. 

From the argument presented above it appears that different formation mecha- 
nisms may apply to different types of filaments. Quiescent filaments and Intermediate 
filaments which mainly fall into the Exterior and Diffuse bipolar region types rely on 
surface effects acting on coronal fields. In contrast for active region filaments a strong 
possibility is the emergence of flux ropes or the formation of flux ropes during emer- 
gence as a result of coronal reconnection. Therefore it is useful to distinguish between 
IF and QF compared to ARF as they may have a different formation mechanism. 



5.6.3 Future Observation of Filament Channels 

The formation, structure, and evolution of solar filaments is an important part of our 
understanding of coronal physics and the behavior of magnetic fields as they are trans- 
ported across the solar surface. New observational results show that surface motions 
acting on non-potential magnetic fields may play an integral part in the formation of 
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large-scale stable filaments while flux rope emergence may play a role for small, unsta- 
ble active region filaments. In addition the large-scale hemispheric pattern of filaments 
may just be the result of surface effects redistributing the helicities of new emerging 
bipoles across the solar surface from low- to high- latitudes. 

A better understanding of the formation of prominences requires multi- wavelength 
observations of prominences situated over a wide range of latitudes, from the active 
region belts up to the polar crowns. It is imperative to determine whether different for- 
mation mechanisms occur at different latitudes on the Sun. To distinguish this, spectral 
lines from Ha to X-rays along with magnetic information are needed to provide full 
coverage of the wavelength ranges associated with the formation and structure of fil- 
aments. Before ordering up this large menu of observations, we must be able to use 
hindsight as to where and when a long-lived filament might form. Therefore, mainte- 
nance of existing synoptic data sets is a vital part of advanced studies of prominence 
formation. 



6 Open Issues 

The present review is focused on four aspects of solar prominences, namely, their mag- 
netic structure, the dynamics of prominence plasma, prominence oscillations, and the 
formation and evolution of filament channels. Other aspects, such as the physical prop- 
erties of prominence plasmas, spectroscopic methods used to measure these properties 
and non-LTE models are discussed in Paper I. In the following we outline several out- 
standing issues, and the observations and modeling required to resolve these issues. 
The ordering reflects the order of topics considered in the review and no attempt is 
made to prioritise them: 

1. Why do different prominences have such different morphology (e.g., horizontal 
vs. vertical threads), and why do prominences/filaments look so different in dif- 
ferent wavelengths? Current models do not readily explain such differences. What 
physical mechanism causes the observed thin thread structures and subsequently 
the gaps between individual threads? To resolve these issues, images with high spa- 
tial resolution of filaments and prominences obtained simultaneously are required 
at different wavelengths, including Ha, He I 10830 A and He II 304 A. These ob- 
servations should track the same filament as it rotates across the solar disk and 
above the limb. Three-dimensional (3D) models of prominence threads should be 
developed (see Paper I). 

2. What is the 3D magnetic structure of different types of prominences? In particular, 
what is the magnetic field orientation in filament barbs, and what is the relation- 
ship between the barbs and the evolving photospheric magnetic fields (e.g., parasitic 
polarities)? What is the effect of the evolving magnetic carpet on prominence mag- 
netic structure? Are non-erupting filaments outside active regions suspended in 
detached flux ropes ? Answering such questions will require spectro-polarimetric 
measurements of the magnetic fields in prominences at multiple heights including 
both measurements in prominences (using the Zeeman and Hanle effects) and in 
the photosphere/chromosphere below filaments. Techniques should be developed 
for modeling the evolution of filament magnetic structure in response to time- 
dependent photospheric boundary conditions, including the motions of discrete 
magnetic elements, flux emergence, and flux cancellation. 
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3. To what extent can vector field extrapolations from the photosphere or chromo- 
sphere be used to deduce magnetic fields within filaments? Answering this ques- 
tion requires photospheric and/or chromospheric vector field measurements, reliable 
NLFFF modeling to extrapolate the magnetic fields to larger heights, and verifica- 
tion that extrapolation results match the observed filament channel structures and 
magnetic fields. 

4. How does the prominence plasma collect and evolve? What physical mechanisms 
are responsible for prominence plasma formation and dynamics? Which mecha- 
nisms best explain the observations of plasma distributions and flows? As different 
magnetic field models exhibit a different field structure (e.g. sheared arcade vs flux 
rope). What effect do these differences have on the distribution and dynamics of 
the plasma produced by the different plasma formation mechanisms. This requires 
simultaneous, high cadence, high spatial resolution observations of many spectral 
lines covering the full range of temperatures present in the solar atmosphere. An 
important goal is to obtain proper motion and Doppler shift information. Tests 
of the proposed mechanisms (§3) requires 3D MHD modeling, including relevant 
energy release, deposition, transport and loss processes. 

5. What is the relationship between the observed signatures of the oscillations (spec- 
tral line shift, width, and intensity variation) and the physical properties of the 
MHD waves (velocity, magnetic field, temperature and density perturbations)? 
How are prominence oscillations excited and damped, and how does this affect the 
plasma energetics? Can prominence seismology be developed as a powerful diagnos- 
tic tool, too for example deduce the distribution and strength of the magnetic field 
in prominences? Solving this problem requires 2D observations of Doppler shifts 
with high spatial and temporal resolution, as well as independent estimates of key 
physical properties (e.g. densities and temperatures) . Interpretation of these obser- 
vations will require multi-dimensional modeling of wave propagation and damping 
in prominences, eventually incorporating wave excitation processes. 

6. How are prominences heated, e.g., by waves (see point 5), shocks or other mecha- 
nisms? The temperatures in prominences are higher than the radiative equilibrium 
temperature of about 5000 K. We need high cadence observations in several spectral 
lines, and accurate modeling of the ionization state of the plasma and of radiative 
transport in optically thick spectral lines (see Paper I). 

7. How are filament channels formed at all latitudes on the Sun? Currently, few obser- 
vations of filament channel formation are available. What is the helicity distribution 
within activity complexes and how does it relate to filament channels? How is this 
helicity dispersed across the solar surface? This requires low cadence (hours) global 
observations, including the far side of the Sun over long periods of time (years): 
photospheric LOS magnetic field, Ha, He I 10830 A, He II 304 A. Techniques should 
be developed to simulate the global-scale, long-term evolution of the coronal mag- 
netic field, incorporating the observed distribution of emerged magnetic flux and 
helicity. 

8. How and at what locations along the PILs do filament channels form within active 
regions and activity complexes? Answering this question requires high-cadence, 
high-resolution observations of active regions and activity complexes in Ha, EUV 
lines, and LOS magnetic field. To elucidate the role of subsurface processes (§5.3 
and §5.6), the emergence of twisted flux ropes should be modeled self-consistently. 

9. What is the magnetic structure of polar crown prominences? Do they have the 
same chirality as mid-latitude prominences in the same hemisphere? If they exhibit 
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a hemispheric pattern what is the cause of it ? This requires spatially resolved, 
synoptic measurements of prominence magnetic fields using the Hanle effect. Such 
observations have not been carried out since the f980's. Do sub-surface processes 
play a role in the formation of polar crown prominences? To answer this question re- 
quires out-of-ecliptic observations of photospheric magnetic fields (LOS component) 
and local helio- seismology to determine subsurface flows near the polar crowns. 
10. What can filaments and prominences tell us about the distribution of magnetic 
helicity across the solar surface and subsequently the solar dynamo? Attacking this 
fundamental issue requires the development of 3D dynamo models that include the 
distribution and transport of magnetic helicity. The results of such dynamo models 
should be compared with observations of the helicity in active regions and filament 
channels on the quiet Sun. 

These issues can be addressed by combining spectral and imaging observations 
from ground- and space-based instruments. Some of the goals may be achieved with 
currently available resources, but others will require new missions. Planners of future 
space missions are encouraged to take the above observational requirements into ac- 
count. 
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